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Abstract 

An  aircraft’s  response  to  control  inputs  varies  widely  throughout  its  flight  envelope.  The 
aircraft  configuration  also  impacts  control  response  through  variations  in  center  of  gravity  and 
moments  of  inertia.  Designing  a  flight  control  system  (FCS)  to  accomodate  the  full  flight  envelope 
and  configuration  set  of  an  aircraft  is  clearly  a  complex  undertaking.  Quantitative  feedback  theory 
(QFT)  is  a  design  tool  which  enables  the  engineer  to  attack  this  task  in  an  efficient  way.  Although 
QFT  is  a  robust  control  design  technique,  it  is  an  interactive  algorithm  allowing  the  engineer  full 
control  over  compensator  order  and  gain.  In  this  research  effort,  a  full  subsonic  flight  envelope  FCS 
is  designed  for  the  VISTA  F-16  aircraft  using  QFT.  Four  aircraft  configurations  are  considered. 
The  strict  control  of  the  compensator  order  and  gain  allowed  by  QFT  facilitates  the  attainment  of 
desired  performance  while  avoiding  physical  saturations.  In  addition,  flying  qualities  are  imbedded 
in  the  longitudinal  design  through  the  use  of  a  control  parameter  which  varies  with  the  aircraft’s 
energy  state.  This  parameter  is  synthesized  to  closely  reflect  the  actual  control  desires  of  the  pilot 
throughout  the  aircraft  flight  envelope.  Linear  simulations  with  realistically  large  control  inputs 
are  used  to  validate  the  design. 
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A  QUANTITATIVE  FEEDBACK  THEORY  PCS  DESIGN  FOR  THE  SUBSONIC 
ENVELOPE  OF  THE  VISTA  F-16  INCLUDING 
CONFIGURATION  VARIATION 

1.  Introduction 

1.1  Research  Justification 

Traditionally,  flight  control  engineers  have  taken  a  conservative,  brute  force  approach  to 
designing  a  full  envelope  flight  control  system  (FCS)  for  an  aircraft.  First,  many  design  points 
within  and  along  the  border  of  the  flight  envelope  of  an  aircraft  are  selected.  These  design  points 
are  usually  airspeed-altitude  combinations  as  shown  by  the  X’s  in  Figure  1.1.  Second,  individual 
compensator  designs  are  accomplished  for  each  of  these  airspeed-altitude  points.  Third,  smooth 
transitions  between  these  compensators  must  be  engineered.  Making  the  transitions  imperceptible 
to  the  pilot  is  very  difficult  and  time-consuming  because  each  airspeed-altitude  design  point  can  be 
approached  from  an  infinite  number  of  initial  conditions.  Obviously,  if  the  number  of  the  airspeed- 
altitude  design  points  can  be  reduced,  thus  reducing  the  number  of  transitions  required,  the  design 
process  can  be  made  more  efficient,  and  the  resulting  FCS  less  complex. 

One  possible  way  to  reduce  the  number  of  necessary  design  points  is  to  apply  a  robust  control 
design  technique  to  the  problem.  A  compensator  synthesized  using  robust  control  principles  should 
certainly  be  able  to  handle  large  parts  of,  if  not  the  whole,  flight  envelope.  Unfortunately,  most 
attempts  at  applying  robust  control  design  algorithms  to  practical,  “real  world”  problems  have  been 
dismal  failures  [7].  This  is  because  although  the  problem  is  well  posed,  the  resulting  compensator 
is  impractical  to  implement.  Either  the  compensator  is  of  too  high  an  order,  or  the  compensator 
gain  is  too  large  to  accommodate  “real  world”  nonlinearities.  Also,  any  sensor  noise  present  is 
accentuated  by  this  high  gain.  The  typical  reason  for  these  poor  results  is  that  the  robust  design 
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I  Airspeed 

Figure  1.1  Airspeed- Altitude  Design  Points  used  in  Traditional  FCS  Design 

is  synthesized  in  the  essentially  noiseless  world  of  the  digital  computer,  and  then  validated  on  the 
digital  computer  through  the  use  of  small  signal,  linear  simulation.  Is  it  any  wonder  that  these 
designs  exhibit  unstable  or  unsatisfactory  behavior  when  removed  from  this  sterilized  environment? 

One  robust  control  design  technique  that  has  the  potential  to  avoid  the  aforementioned  pit- 
falls  is  quantitative  feedback  theory  (QFT).  QFT  is  a  mathematically  rigorous,  frequency  domain, 
robust  control  design  technique.  Although  a  QFT  design  effort  could  very  easily  result  in  a  com¬ 
pensator  of  high  order  and  of  high  gain,  it  does  give  the  designer  complete  control  over  the  gain  and 
order  of  the  compensator;  hence,  QFT  is  not  constrained  to  produce  an  impractical  compensator. 
In  addition,  if  a  decision  is  made  to  decrease  or  limit  the  order  or  gain  of  a  compensator,  the 
performance  tradeoffs  due  to  this  action  can  be  clearly  seen  by  the  designer. 

Due  to  the  complex  mathematics  involved,  the  QFT  design  process  can  be  very  tedious  and 
time-consuming  if  applied  without  the  help  of  a  capable  computer-aided  design  (CAD)  tool.  As  a 
result,  most  QFT  design  efforts  prior  to  the  development  of  a  QFTCAD  package  by  Sating  in  1992 
were  constrained  to  be  of  very  limited  scope  [9]. 
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The  only  full  subsonic  flight  envelope  design  effort  using  QFT  has  been  that  of  Reynolds  in 
1993  [8].  Reynolds  used  QFT  to  design  an  FCS  for  the  subsonic  flight  envelope  of  the  VISTA  F-16. 
Although  Reynolds  included  all  variations  in  the  aircraft’s  response  to  control  surface  deflections 
due  to  airspeed  and  altitude  changes,  he  did  not  include  any  variations  caused  by  changes  in  the 
aircraft  center  of  gravity  (CG).  Center  of  gravity  changes  are  a  result  of  either  fuel  consumption 
or  the  addition  or  expenditure  of  external  stores.  Reynold’s  recommendations  were  to  extend  his 
research  to  include  CG  variations  along  with  those  caused  by  airspeed  and  altitude  changes  [8], 

In  summary,  although  excellent  FCS’s  have  been  designed  for  aircraft  using  traditional  design 
methods,  the  synthesis  of  those  FCS’s  has  been  a  costly,  time-consuming  endeavor.  Thus,  limiting 
robustness  in  FCS  design  results  in  a  convoluted,  complex,  full  envelope  design.  QFT  offers  the  hope 
of  incorporating  enough  robustness  to  simplify  the  design  process  and  the  resulting  FCS,  but  not 
so  much  robustness  that  the  resulting  FCS  is  impractical  to  implement  due  to  violation  of  physical 
limitations  imposed  by  the  “real  world”  (i.e.,  actuator  saturation  or  sensor  noise  amplification). 

l.S  Research  Objectives 

The  objectives  of  this  research  effort  are; 

1.  Design  a  practiced  FCS  for  the  subsonic  flight  envelope  of  the  VISTA  F-16  using  QFT. 

(a)  Include  parameter  uncertainty  due  to  variations  in  the  aircraft  CG. 

(b)  Keep  compensator  gain  scheduling  to  an  absolute  minimum. 

(c)  Achieve  performance  comparable  to  the  current  F-16  FCS. 

(d)  Achieve  Level  1  flying  qualities. 

(e)  Critically  validate  the  resulting  design  with  full  consideration  of  “real  world”  factors 
such  as  control  surface  deflection  limits  and  actuator  rate  limits. 
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2.  Further  QFT’s  standing  as  a  viable  robust  control  design  method  that  can  be  applied  to 


practical  problems  and  produce  implementable  results. 
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II.  Research  Methodology 


The  primary  tool  used  to  accomplish  this  flight  control  system  design  is  QFT.  However, 
due  to  the  necessity  to  keep  in  touch  with  the  “real  world”,  an  integrated  design  process  which 
includes  time  domain  simulation  is  utilized.  Although  such  a  process  involves  some  iteration,  the 
flexibility  and  insight  provided  by  QFT  keep  it  to  a  minimum.  In  other  words,  once  a  violation 
of  a  physical  limitation  is  identified,  the  designer  merely  makes  an  insightfull  adjustment  to  the 
present  design  rather  than  repeating  another  complete  design  iteration.  Although  unacceptable  in 
mathematics,  such  trial  and  error  is  an  accepted  engineering  practice.  The  parameter  uncertainty 
for  the  QFT  design  is  quantified  by  linear  time-invariant  perturbation  aircraft  models  collected  for 
many  airspeed-altitude  data  points  throughout  and  along  the  edges  of  the  VISTA  F-16  subsonic 
flight  envelope.  These  aircraft  models  are  generated  for  several  different  aircraft  configurations. 
Generation  of  the  aircraft  models  is  discussed  in  Chapter  III.  Frequency  domain  validation  is 
performed  by  the  QFTCAD,  and  time  domain  validation  is  accomplished  by  examining  the  results 
of  Matlab  linear  time  simulations  of  the  system’s  response  to  realistically  large  input  signals. 

2.1  QFT  Description 

A  thorough  discussion  of  the  QFT  design  process  is  presented  in  the  references  and  is  not 
reproduced  here  [5].  However,  the  fundamentals  of  the  technique  as  they  apply  to  this  research  are 
addressed. 

QFT  is  a  two  degree  of  freedom  design  technique  utilizing  the  fundamental  multiple-input 
single-output  (MISO)  structure  illustrated  in  Figure  2.1  [5,6].  The  two  degrees  of  freedom  available 
to  the  designer  are  the  compensator  G  and  the  prefilter  F.  The  characteristics  of  the  uncompensated 
system  are  contained  in  the  plant  P. 

2.1.1  Quantifying  Parameter  Uncertainty.  Before  beginning  a  QFT  design,  the  un¬ 
certainty  to  be  handled  by  the  design  must  be  quantified.  In  the  case  of  an  FCS  design,  this 
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Disturbance  D 


Figure  2.1  QFT  MISO  Design  Structure 


uncertainty  is  the  variation  in  aircraft  response  to  control  surface  deflections  caused  primarily  by 
changes  in  airspeed,  altitude,  and  CG  position.  This  uncertainty  is  embodied  in  the  aircraft’s 
stability  derivatives.  Therefore,  to  quantify  the  uncertainty,  aircraft  stability  derivatives  for  many 
airspeed-altitude-CG  combinations  spanning  the  flight  envelope  and  aircraft  configurations  of  in¬ 
terest  are  collected.  These  stability  derivatives  are  assembled  into  transfer  functions  representing 
the  aircraft  response  to  the  deflection  of  its  control  surfaces.  These  input-output  transfer  functions 
are  referred  to  as  the  aircraft  plants. 

2.1.2  Frequency  Templates.  Frequency  templates  are  used  in  QFT  to  display  parameter 
uncertainty.  A  frequency  template  is  generated  by  plotting  the  open-loop  frequency  response  of 
each  aircraft  plant  at  a  particular  frequency.  The  result  is  a  collection  of  points  (one  for  each 
plant)  on  a  two-dimensional  graph  of  relative  angle  versus  relative  magnitude  as  illustrated  by  the 
X’s  in  Figure  2.2.  Then,  as  shown  in  the  figure,  a  polygon  is  drawn  by  connecting  the  outermost 
frequency  response  points  so  that  all  other  points  are  contained  within  the  polygon.  This  polygon  is 
the  template  which  corresponds  to  a  given  frequency  and  displays  the  uncertainty,  in  the  frequency 
domain,  embodied  by  all  of  the  aircraft  plants  at  that  particular  frequency.  A  family  of  templates 
is  generated  for  a  representative  set  of  frequencies  in  the  range  of  interest  for  a  particular  design. 
The  resulting  family  of  templates  quantifies  the  parameter  uncertainty  of  the  system  across  the  full 
design  frequency  spectrum. 

The  overall  size  of  the  templates  indicates  to  the  designer  how  much  uncertainty  the  compen¬ 
sator  will  have  to  handle  and  whether  or  not  a  single  compensator  will  be  able  to  do  the  job.  In 
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Figure  2.2  QFT  Frequency  Template 

addition,  if  it  is  determined  that  the  compensator  design  must  be  split  into  more  than  one  part  (i.e., 
full  envelope  robustness  cannot  be  achieved),  the  templates  provide  valuable  insight  as  to  where  in 
the  parameter  space  the  design  should  be  divided. 

2.1.3  Definition  of  Specifications.  In  a  QFT  design,  it  is  necessary  to  define  frequency 
domain  specifications  for  tracking  response,  stability,  and  external  disturbance  rejection. 

Design  tracking  specifications  are  typically  defined  in  the  time  domain  and  must  be  trans¬ 
formed  into  the  frequency  domain  in  order  to  be  used  in  QFT.  This  transformation  is  accomplished 
by  modelling  the  acceptable  range  of  desired  time  domain  responses  of  the  system  by  two  transfer 
functions  (e.g.,  one  representing  the  most  underdamped  response  allowed,  and  the  other  represent¬ 
ing  the  most  overdamped  response  allowed).  The  frequency  response  of  these  models  is  plotted  to 
form  upper  and  lower  tracking  bounds,  Tru  and  Trr,  as  depicted  in  Figure  2.3.  The  difference 
between  Tru  and  Trl  at  a  particular  frequency  is  referred  to  as  6r  and  represents  the  acceptable 
variation  between  the  frequency  responses  of  all  of  the  plants  at  that  frequency.  In  other  words,  6r 
defines  the  frequency  domain  robustness  required  of  the  design  at  each  frequency.  To  facilitate  the 
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Figure  2.3  QFT  'Racking  Specification  Bounds 


design  of  the  compensator,  Tru  and  Trc  should  be  synthesized  so  that  Sr  monotonically  increases 
above  the  0  decibel  (dB)  crossing  frequency  ofTfli;  [5], 

Stability  specifications  are  typically  expressed  in  frequency  domain  terms  of  open-loop  phase 
and/or  gain  margins  and  can  be  used  directly  in  QFT. 

External  disturbance  rejection  specifications  for  a  disturbance  input  to  the  plant,  as  shown  by 
D  in  Figure  2.1,  are  typically  modelled  by  a  constant  value  (e.g.,  -20  dB)  for  all  frequencies  as  shown 
by  Tdu  in  Figure  2.4.  Since  disturbances  to  the  plant  input  are  the  only  disturbance  of  interest  in 
this  FCS  design,  disturbances  to  the  plant  output  are  not  discussed;  however,  disturbances  to  the 
output  can  be  accounted  for  in  a  QFT  design  if  required  [5] . 


2.1.4  QFT  Design  Nichols  Chart  Boundaries.  The  uncertainty  embodied  by  the  fre¬ 
quency  templates  is  used  along  with  the  desired  performance  dictated  by  the  tracking,  stability, 
and  disturbance  rejection  models  to  generate  QFT  design  boundaries  on  a  Nichols  chart.  These 
design  boundaries  are  used  to  guide  the  synthesis  of  a  compensator  to  achieve  the  desired  results. 


2-4 


Frequency  (radians/sec) 


♦ 


Figure  2.4  QFT  External  Disturbance  Rejection  Specification  Bound 


Prior  to  generating  these  design  boundaries,  the  designer  selects  a  nominal  plant.  The  system 
open-loop  transfer  function  (OLTF)  generated  using  this  nominal  plant  is  the  only  one  used  in  the 
design  of  the  compensator.  The  nominal  plant  is  typically  a  plant  located  at  a  particular  position 
along  the  border  of  the  frequency  templates.  This  selection  is  heavily  dependent  on  the  specific 
design  requirements,  but  some  general  guidelines  are  presented  in  the  reference  [5]. 

Once  a  nominal  plant  is  selected,  Nichols  chart  tracking,  external  disturbance  rejection,  and 
stability  design  boundaries  are  generated  for  each  one  of  the  template  frequencies.  Referencing 
Figure  2.1,  the  nominal  OLTF  from  the  output  of  the  prefilter  F  to  the  system  output  y  is 


Lo  =  GP„ 


(2.1) 


where  Po  is  the  nominal  plant.  The  design  boundaries  are  generated  to  direct  the  shaping  of  this 
nominal  loop,  Lg. 

The  previously  computed  frequency  templates  are  used  to  synthesize  the  Nichols  chart  bound¬ 
aries,  which,  if  satisfied  by  the  nominal  loop,  guarantee  that  all  of  the  open-loop  functions  £-,•  =  GP,- 
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Figure  2.5  QFT  Nichols  Chart  Design  Boundaries 


meet  the  tracking,  stability,  and  external  disturbance  rejection  specifications.  The  procedure  used 
to  generate  these  boundaries  is  thoroughly  discussed  in  the  reference  [5].  An  example  of  these 
Nichols  chart  boundaries  is  shown  in  Figure  2.5,  where  the  stability  boundaries  are  the  closed 
boundaries  encircling  the  {0  dB,  -180  degree}  point,  and  the  horizontal  boundaries  are  a  composite 
boundary  formed  by  displaying  the  most  restrictive  (i.e.,  the  highest  on  the  Nichols  chart)  of  the 
tracking  and  disturbance  rejection  boundaries  at  each  phase  angle  [1].  Note  that  there  is  a  stability 
and  composite  boundary  for  each  one  of  the  template  frequencies. 

2.1.5  QFT  Compensator  Design.  The  compensator  G  of  Figure  2.1  is  now  synthesized 
by  plotting  Lg  =  GPg  on  the  Nichols  chart  containing  the  design  boundaries.  Poles  and  zeros  are 
added  to  G  and  its  gain  adjusted  until  the  boundaries  are  satisfied  by  the  nominal  loop.  If  the 
uncompensated  OLTF  Pg  is  not  type  1  (i.e.,  it  lacks  pure  integral  action),  the  first  step  taken  is 
to  give  G  a  pole  at  the  origin  to  guarantee  zero  steady-state  tracking  error  of  a  step  input  and  to 
provide  disturbance  rejection  [5]. 
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To  satisfy  a  tracking  or  disturbance  boundary,  the  magnitude-phase  point  of  corresponding 
to  the  boundary  frequency  must  lie  on  or  above  the  boundary.  Also,  to  satisfy  a  stability  boundary 
the  magnitude-phase  point  of  Lg  corresponding  to  the  boundary  frequency  must  lie  outside  the 
boundary.  A  typical  loop  shaping  Nichols  chart  is  depicted  in  Figure  2.6  [1].  The  nominal  loop  is 
the  line  which  starts  at  the  top  of  the  figure  at  -180  degrees  of  phase  and  then  curves  right  to  avoid 
the  stability  boundaries  and  finally  curves  left  once  below  the  stability  boundaries. 


Satisfying  the  tracking  boundary  at  a  particular  frequency  ensures  that  the  required  robust¬ 
ness  is  achieved  by  all  of  the  closed-loop  frequency  responses  at  that  frequency.  In  other  words, 
the  difference  between  the  smallest  magnitude  and  largest  magnitude  frequency  response  of  all  the 
closed-loop  functions 


GPj 
1  +  GPi 


(2.2) 


is  less  than  or  equal  to  the  6jt  (defined  in  Figure  2.3)  that  corresponds  to  that  tracking  boundary 
frequency. 

The  synthesis  of  a  compensator  that  ensures  the  nominal  loop  satisfies  all  of  the  tracking 
boundaries  results  in  closed-loop  frequency  responses  like  those  shown  in  Figure  2.7.  As  is  seen  in 
the  figure,  all  of  the  closed-loop  frequency  responses  are  within  Sn  of  one  another  at  all  frequencies, 
but  they  do  not  yet  fall  between  Tjiu  and  Trl  ;  hence,  although  the  required  robustness  has  been 
achieved,  the  performance  of  the  system  does  not  yet  meet  specifications. 

Satisfying  the  disturbance  boundary  at  a  particular  frequency  ensures  that  the  magnitude 
of  all  of  the  closed-loop  frequency  responses  at  that  frequency  from  the  disturbance  input  D  of 
Figure  2.1  are  less  than  Tou  of  Figure  2.4.  It  follows  that  if  a  compensator  is  synthesized  which 
ensures  Lg  satisfies  all  of  the  disturbance  rejection  boundaries,  then  the  frequency  responses  of  all 
the  disturbance  input  closed-loop  functions 


Td.= 


Pi 

1  +  GPi 


(2.3) 
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Figure  2.6  QFT  Nominal  Loop  Shaping  Nichols  Chart 


Figure  2.7  Tracking  Input  Closed-Loop  Frequency  Responses  after  Compensator  Design 
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Figure  2.8  Disturbance  Input  Closed-Loop  Frequency  Responses  after  Compensator  Design 

lie  below  Tdu  for  all  frequencies  as  shown  in  Figure  2.8. 

If  the  composite  boundaries  are  satisfied,  then  the  results  displayed  in  Figures  2.7  and  2.8  are 
achieved  simultaneously. 

The  necessary  stability  of  the  system  can  only  be  achieved  by  shaping  Lo  to  remain  outside 
the  closed  stability  boundaries  at  all  frequencies. 

This  is  the  point  in  the  QFT  design  process  at  which  the  designer  can  evaluate  the  gain 
and  order  of  G.  If  G  needs  to  be  reduced  in  order  or  gain,  the  designer  can  graphically  see  the 
design  tradeoffs  involved.  The  tradeoffs  that  can  be  readily  evaluated  on  the  Nichols  chart  include 
degredation  of  tracking  response  or  disturbance  rejection,  decrease  of  stability  (phase  margin  angle 
or  gain  margin  dB),  and  the  increase  in  phase  or  gain  margin  frequencies.  Therefore,  the  designer 
can  not  only  intelligently  reduce  the  order  or  gain  of  G,  but  is  immediately  aware  of  the  consequences 
of  doing  so. 


2.1.6  QFT  Prefilter  Design.  The  purpose  of  the  prefilter  F  of  Figure  2.1  is  to  shift  the 
robust  frequency  responses  of  Figure  2.7  so  that  they  fall  between  Tru  and  Trl-  Therefore,  this 


Figure  2.9  Successful  QFT  Design  Closed- Loop  Frequency  Responses 


final  step  in  the  QFT  design  process  results  in  a  design  that  produces  the  desired  performance. 
Figure  2.9  illustrates  the  frequency  responses  of  all  the  closed-loop  functions 


FGPi 

•  "  1  -1-  GPi 


(2.4) 


resulting  from  a  successful  QFT  design. 

2.1.7  Multiple-Input  Multiple- Output  QFT.  All  of  the  preceding  discussion  of  the  QFT 
design  process  has  been  directed  towards  MISO  systems,  but  QFT  can  also  be  applied  to  multiple- 
input  multiple-output  (MIMO)  control  systems  like  the  one  in  Figure  2.10.  This  is  possible  because 
any  MIMO  system  can  be  represented  equivalently  as  rrp  MISO  unity  feedback  loops,  where  m  is 
the  dimension  of  a  square  MIMO  system.  The  full  proof  of  this  equivalence  can  be  found  elsewhere 
and  only  the  definition  of  the  equivalent  MISO  loops  will  be  presented  here  [5,6]. 

The  structure  of  the  equivalent  MISO  loops  of  a  2  X  2  MIMO  control  system  is  shown  in 
Figure  2.11.  The  first  subscript  used  in  the  figure  refers  to  the  output  and  the  second  subscript 
refers  to  the  input.  As  an  example,  j/12  refers  to  the  part  of  the  first  output  of  the  MIMO  system 
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Figure  2.10  2X2  MIMO  Control  System 


MISO  Loop  {2,11  MISO  Loop  (2^) 


Figure  2.11  Equivalent  QFT  MISO  Loops  of  a  2  X  2  MIMO  Control  System 


that  is  a  result  of  an  input  to  the  second  MIMO  input  port.  The  prefilter  and  compensator  are 
typically  designed  as  diagonal  and  hence  have  the  following  structure. 


F 


G 


h  0 
0  /2 

ffi  0 
0  52 


(2.5) 


(2.6) 


The  Cij’s  in  Figure  2.11  represent  the  internal  cross-coupling  interaction  of  all  the  other  loops 
with  the  loop  of  interest.  The  presence  of  the  Cj/s  in  the  design  structure  makes  this  MISO  structure 
fully  representative  of  the  original  MIMO  system. 
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The  g,j’s  of  Figure  2.11  are  the  elements  of  the  equivalent  MISO  structure  that  embody  the 
characteristics  of  the  plant  and  are  defined  as  follows.  Let  the  inverse  of  the  original  plant  be 
defined  as 


Pll* 

Pl2 

P21* 

P22' 

The  qij ’s  are  then  defined  as  the  elements  of  the  matrix 


1 

1 

Pll* 

PI2* 

1 

1  * 

Pai* 

P32 

(2.7) 


(2.8) 


Once  the  Q  matrix  has  been  calculated,  each  of  the  four  loops  can  be  designed  using  the  MISO 
procedures  outlined  in  the  previous  sections. 


2.1.8  QFTCAD.  Fortunately,  the  QFTCAD  package  developed  by  Sating  automates  the 
design  tasks  described  in  the  preceding  sections  amd  is  utilized  heavily  during  this  design  effort  [1]. 
Without  this  package,  a  QFT  design  of  this  magnitude  would  be  prohibitively  tedious  and  time- 
consuming. 


2.2  Time  Simulation  Description 

All  time  simulations  in  this  research  are  performed  using  Matlab’s  ‘rk45’  algorithm  [11].  This 
algorithm  performs  Runge-Kutta  fifth  order  integration  using  fourth  order  step-size  control.  The 
tolerance  of  the  integration  is  set  at  10“®  and  the  minimum  step  size  is  set  to  10“*. 

Also,  the  time  simulations  are  conducted  for  a  period  of  only  five  seconds  for  two  reasons. 
First,  when  controlling  an  aircraft,  pilots  make  control  inputs  more  frequently  than  once  every  five 
seconds.  To  consider  simulation  for  longer  than  a  five  second  period  would  not  make  sense  since 
the  pilot  will  most  likely  have  made  several  control  inputs  during  that  period.  In  other  words,  the 
low  end  of  a  pilot’s  control  bandwidth  is  greater  than  0.2  radians  per  second  (rps).  Second,  the 
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linear  perturbation  aircraft  models  used  in  this  design  do  not  remain  valid  for  simulations  extended 
beyond  five  seconds. 

2.3  Specific  Design  Guidelines 

During  the  QFT  longitudinal  and  lateral  channel  compensator  synthesis,  several  guidelines  are 
followed.  First,  the  primary  frequency  range  of  interest  for  meeting  tracking  response  specifications 
is  0.5  to  3.5  rps,  which  corresponds  to  a  typical  pilot’s  control  bandwidth  [7]. 

Second,  the  set  of  template  frequencies  used  during  the  QFT  portion  of  the  design  is 

w,-  =  {0.05,  0.075,  0.1,  0.25,  0.5,  0.75,  1,  3,  5,  7,  10,  15,  20,  25,  30,  50}  rps.  (2.9) 

Third,  the  compensators  are  designed  so  that  the  bending  modes  of  the  aircraft  are  not 
excited.  These  bending  modes  are  isolated  to  frequencies  above  30  rps;  hence,  to  prevent  their 
excitation,  the  phase  margin  frequency  of  all  the  compensated  OLTF’s  must  be  kept  at  or  below 
30  rps  [7],  To  satisfy  this  requirement,  the  nominal  plant  for  this  design  is  chosen  as  the  plant 
located  at  the  top  of  the  30  rps  frequency  template.  In  this  case,  if  the  nominal  OLTF  has  an 
of  30  rps,  it  is  guaranteed  that  all  other  plant  OLTF’s  have  an  of  less  than  30  rps;  therefore, 
none  of  the  plant  responses  interact  with  the  bending  modes  of  the  aircraft.  This  choice  of  nominal 
plant  is  the  only  one  that  allows  the  designer  to  easily  limit  the  of  all  plant  OLTF’s. 

Fourth,  restrictions  are  placed  on  the  compensators  themselves.  The  resulting  compensators 
must  be  of  third-order  or  less,  and  no  poles  of  the  compensators  are  to  be  faster  than  60  rps.  The 
first  of  these  restrictions  is  to  keep  the  compensator  dimension  to  an  implementable  level.  The 
second  restriction  is  imposed  in  anticipation  of  future  digitization  of  the  design.  The  detrimental 
alictsing  effects  present  when  digitizing  at  the  typical  sampling  rates  used  in  flight  control  (e.g.,  50 
Hz)  are  avoided  if  the  compensator  poles  are  selected  to  comply  with  the  60  rps  restriction. 
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Control  Surface  i 

Deflection  Limit 

horizontal  tail 

±20  degrees 

flaperons 

±20  degrees 

rudder 

±30  degrees 

Table  2.1  Control  Surface  Deflection  Limits 


Fifth,  the  physical  limitations  of  actual  aircraft  hardware  that  restrict  the  FCS  design  are 
deflned.  One  piece  of  hardware  considered  is  the  control  surface  hydro-mechanical  actuator.  The 
actuator  has  a  60  degree  per  second  rate  limit.  If  a  rate  higher  than  this  is  demanded,  the  actuator  is 
saturated.  The  aircraft  control  surfaces  also  impose  restrictions  on  the  FCS  design.  The  deflection 
limit  used  in  this  design  for  each  of  the  aircraft  control  surfaces  is  listed  in  Table  2.3.  If  the 
FCS  demands  more  deflection  than  is  physically  available,  control  surface  saturation  occurs.  Both 
actuator  rate  and  control  surface  deflection  saturation  may  significantly  degrade  performance  and 
reduce  FCS  stability,  and  are  therefore  avoided  in  this  design. 

Finally,  the  prefilters  chosen  must  not  restrict  the  control  bandwidth  of  the  pilot.  The  typical 
bandwidth  of  a  pilot  is  0.5  to  3.5  rps  [7];  Consequently,  the  prefilters  of  this  design  must  have 
bandwidths  of  no  less  than  3.5  rps. 
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III.  FCS  Design  Aircraft  Model 


To  quantify  the  parameter  uncertainty  involved  in  a  subsonic  flight  envelope  FCS  design 
including  changes  in  the  aircraft  CG,  linear  time-invariant  (LTI)  aircraft  models  are  generated  that 
span  the  airspeed-altitude-CG  parameter  space.  These  aircraft  models  are  used  as  the  plants  in 
the  ensuing  QFT  design. 

3.1  LTI  Aircraft  Model  Generation 

The  tool  chosen  to  generate  the  required  LTI  aircraft  models  is  the  Simulation/Rapid- 
Prototyping  Facility  (SRF)  which  is  available  at  the  Wright-Patterson  Air  Force  Base  Flight  Dynam¬ 
ics  Laboratory  [4].  Once  an  altitude,  airspeed,  and  configuration  are  input  to  the  SRF,  nonlinear 
equations  of  motion  are  linearized  about  this  operating  point  and  stability  derivatives  in  the  form 
of  state  space  A  and  B  matrices  are  generated.  The  arrangement  of  the  stability  derivatives  in  the 
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state  space  equation  x  —  Ax  +  Bu  is  shown  in  Equation  (3.1). 
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(3.1) 


The  units  of  each  of  the  states  listed  in  Equation  (3.1)  are  listed  in  Table  3.1.  It  is  evident  from 
this  state  space  representation  that  the  longitudinal  and  lateral  channels  are  considered  completely 
decoupled  or  independent  of  one  another.  This  is  consistent  with  the  steady  level  trim  conditions 
at  which  the  plants  are  generated. 
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State 

Units 

e 

deg 

u 

ft /sec 

a 

deg 

9 

deg/sec 

<i> 

deg 

0 

deg 

p 

deg /sec 

r 

deg/sec 

Table  3.1  Units  of  Aircraft  Model  States 


3.S  Longitudinal  Aircraft  Model 

The  state  space  longitudinal  aircraft  model  extracted  from  Equation  (3.1)  is 
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The  flap  function  of  the  F-16  flaperons  is  selected  only  during  the  takeoff  and  landing  phases  of 
flight;  hence,  since  the  takeoff  and  landing  phases  of  flight  are  not  considered  in  this  PCS  design, 
the  Sfiap  control  input  has  been  eliminated. 

Pilot  control  inputs  are  typically  spaced  no  further  apart  than  five  seconds;  therefore,  five 
seconds  is  accepted  as  an  appropriate  time  period  of  interest  in  flight  control  system  design.  Due 
to  the  dominance  of  the  short  period  mode  in  the  five  seconds  following  a  control  input,  the  short 
period  approximation  is  used  in  the  longitudinal  channel  design  process.  The  resulting  short  period 
approximation  state  space  model  is  shown  in  Equation  (3.3). 
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In  Equation  (3.3),  the  control  input  SeUv  refers  to  the  symmetric  component  of  aircraft  horizontal 
tail  deflections. 

The  corresponding  state  space  output  equation  y  =  Cx+Du  for  this  short  period  longitudinal 
system  is, 


One  final  action  taken  is  to  change  the  sign  of  the  SRF  generated  B  matrix  of  Equation  (3.3). 
This  action  establishes  the  convention  of  positive  elevator  deflection  resulting  in  positive  changes 
in  q  and  a. 

In  summary,  since  the  longitudinal  channel  is  fully  decoupled  from  the  lateral  channel  (see 
Equation  (3.1))  and  the  flaps  are  eliminated  as  a  control  input,  the  resulting  longitudinal  model 
represents  a  single-input  single-output  (SISO)  system  void  of  any  cross-coupling  effects  from  the 
lateral  channel. 

3.3  Lateral  Aircraft  Model 

Unlike  the  longitudinal  channel  where  an  approximation  (short  period)  is  made,  the  full  lateral 
state  space  model  including  the  spiral,  dutch  roll,  and  roll  modes  is  used  in  the  lateral  channel 


3-4 


design.  The  lateral  state  space  model  extracted  from  Equation  (3.1)  is  listed  in  Equation  (3.5). 
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In  Equation  (3.5)  the  three  control  inputs  Sdjtaii  >  ^aii  .  and  drud  correspond  to  differential  deflection 
of  the  horizontal  tails,  aileron  deflection  of  the  flaperons,  and  rudder  deflection,  respectively. 


The  corresponding  state  space  output  equation  y  —  Cx  Du  for  the  lateral  system  is. 
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A  final  action  taken  is  to  change  the  sign  of  the  SRF  generated  B  matrix  of  Equation  (3.5). 
This  action  establishes  the  convention  of  positive  aileron  and  differential  tail  deflections  resulting  in 
positive  changes  in  p  and  the  convention  of  positive  rudder  deflection  resulting  in  positive  changes 
in  r  and  negative  changes  in  /?. 

The  resulting  lateral  model  represents  a  3  X  2  MIMO  system  void  of  any  cross-coupling  effects 
from  the  longitudinal  channel. 
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3-4  Parameter  Space  Boundaries 


The  subsonic  flight  conditions  spanned  by  this  design  include  all  flying  airspeeds  between  and 
including  the  altitudes  of  1,000  and  50,000  feet.  Variation  in  aircraft  CG  position  is  accomplished  by 
generating  models  for  the  preceding  flight  conditions  for  each  of  four  different  aircraft  configurations 
available  in  the  SRF  [4] .  These  configurations  are: 

1.  Clean  (no  external  fuel  tanks) 

2.  Centerline  fuel  tank 

3.  Two  wing  tanks 

4.  Three  external  fuel  tanks  (centerline  and  wings) 

In  Reynold’s  thesis,  the  first  of  these  configurations  was  the  only  one  considered. 

Reynold’s  procedure  for  determining  the  low  airspeed  (left)  side  of  the  flight  envelope  by 
examining  the  behavior  of  the  short  period  and  phugoid  roots  of  the  aircraft  longitudinal  model  of 
Equation  (3.2)  is  repeated  in  this  thesis  [8].  The  addition  of  external  fuel  tanks  does  have  an  effect 
on  the  airspeed  below  which  the  SRF  will  initialize  and  on  the  airspeed  at  which  the  short  period 
and  phugoid  mode  poles  become  well  behaved.  Hence,  the  left  edge  of  the  flight  envelope  is  slightly 
different  for  each  configuration. 

Another  consequence  of  the  external  fuel  tanks  is  the  introduction  of  transonic  effects  at  a 
lower  airspeed  than  for  a  clean  configuration.  Therefore,  the  behavior  of  the  short  period  and 
phugoid  mode  poles  between  0.8  and  0.9  Mach  for  the  configurations  with  fuel  tanks  is  observed. 
The  highest  Mach  at  which  the  modes  become  well  behaved  is  used  as  the  high  airspeed  (right) 
side  of  the  subsonic  flight  envelope. 
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3.5  Validation  of  Parameter  Space  Data  Points 


One  important  issue  to  be  settled  is  whether  or  not  any  flight  condition/configuration  data 
points  exist  in  the  parameter  space  that  contain  plant  variation  not  already  accounted  for  by  the 
data  points  already  selected.  In  other  words,  it  needs  to  be  validated  that  every  bit  of  the  aircraft 
plant  variation  for  the  design  parameter  space  is  contained  in  the  generated  LTI  models.  This  can 
be  determined  by  examining  QFT  frequency  templates. 

First,  all  of  the  short  period  longitudinal  transfer  functions  from  to  pitch  rate  q  are 
loaded  into  the  QFTCAD.  The  details  of  the  procedure  for  generating  transfer  functions  using 
Matlab  and  then  formatting  them  for  input  into  the  QFTCAD  are  included  in  Appendix  A.  This 
generation  procedure  is  used  throughout  this  research  in  order  to  handle  a  large  number  of  plants 
efficiently. 

Second,  the  plants  that  form  the  perimeter  of  each  frequency  template  are  noted  from  an 
expanded  template  view  like  that  shown  in  Figure  3.1  [1].  In  all,  66  plants  are  labelled  as  perimeter 
plants.  Next,  the  flight  conditions  corresponding  to  these  perimeter  plants  are  referenced.  If  all 
of  the  perimeter  plants  correspond  to  flight  condition/configuration  data  points  on  the  edge  of 
the  flight  envelope,  then  the  templates  can  only  be  expanded  by  gathering  data  for  invalid  flight 
conditions  outside  the  actual  VISTA  F-16  flight  envelope. 

After  referencing  the  flight  conditions  used  by  the  SRF  to  generate  the  perimeter  plants,  it  is 
determined  that  all  of  these  plants  are  “edge  of  the  envelope”  plants  except  for  a  handful  between 
the  altitudes  of  10,000  and  30,000  feet  and  between  the  airspeeds  of  0.5  and  0.7  Mach.  Consequently, 
additional  plants  are  generated  by  the  SRF  using  altitudes  of  10,000,  20,000,  and  30,000  feet  and 
airspeeds  at  0.05  Mach  increments  between  0.4  and  0.8  Mach.  The  short  period  approximation 
of  the  resulting  plants  is  added  to  the  66  perimeter  plants  and  input  into  the  QFTCAD  to  see  if 
any  increase  in  template  size  is  achieved.  Figures  3.2  and  3.3  depict  the  templates  for  the  original 
perimeter  plants  and  those  for  the  original  plus  the  added  plants,  respectively.  The  difference 
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Figure  3.1  Expanded  View  of  Frequency  Template 


between  the  two  sets  of  templates  is  imperceptible  in  these  figures.  An  expanded  template  view 
reveals  that  a  slight  increase  in  template  size  (lower  left  region)  is  achieved  by  adding  the  extra 
plants.  This  expanded  view  is  illustrated  for  a  frequency  of  1  rps  in  Figure  3,4.  The  numbers  above 
66  in  the  figure  represent  the  added  plants.  Also,  the  example  expanded  template  in  Figure  3,1  is 
the  template  for  only  the  original  66  perimeter  plants  The  increase  in  template  size  can  be  seen 
when  it  is  compared  to  the  template  of  Figure  3.4.  Since  the  additional  plants  do  increase  the  size 
of  the  frequency  templates,  they  are  retained  as  parameter  space  data  points. 

At  the  completion  of  this  QFT  frequency  template  parameter  space  validation,  it  is  stated  with 
some  certainty  that  the  frequency  templates  do  encompass  all  of  the  plant  uncertainty  introduced 
by  variations  in  flight  condition  and  configuration  throughout  the  full  subsonic  flight  envelope  of 
the  VISTA  F-16. 

A  complete  list  of  the  flight  condition/configuration  data  points  used  in  the  design  is  contained 
in  Appendix  B.  In  addition,  a  tabulation  of  SRF-generated  stability  derivatives  and  other  pertinent 
model  data  is  found  in  the  supplement  to  this  thesis. 
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Figure  3.2  QFT  Frequency  Templates  of  Original  Perimeter  Plants 


Figure  3.3  QFT  Frequency  Templates  of  Original  and  Added  Validation  Plants 
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Figure  3.4  Expanded  View  of  Frequency  Template  for  w  =  1  rps 
3.6  Control  Surface  Actuator  Model 

The  fourth-order  actuator  model  used  in  Reynold’s  thesis  is  also  chosen  for  this  design  and 
has  the  transfer  function  given  in  Equation  (3.7). 


Sccntroijs)  _ (20.2)(71.4)^(144.8) _ 

(s)  (s  -f  20.2)(s  -|-  144.8)[s2  +  2(0.736)(71 .4)s  +  (71.4)2] 


The  “control”  subscript  in  Equation  (3.7)  refers  to  any  of  the  control  surfaces  of  the  aircraft.  As 
discussed  in  Reynold’s  thesis,  a  fourth-order  model  is  chosen  instead  of  a  first-order  model  due  to 
its  superior  modelling  of  the  actual  actuator  frequency  response  phase  at  higher  frequencies  [8]. 
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IV.  Longitudinal  Channel  Design 

An  integrated  approach  is  used  to  accomplish  the  longitudinal  channel  FCS  design.  Quanti¬ 
tative  feedback  theory  is  used  to  design  a  compensator  that  meets  frequency  domain  specifications 
and  then  time  domain  simulation  is  used  to  determine  if  any  control  surface  deflection  or  actuator 
rate  limits  are  exceeded.  If  either  of  these  physical  limitations  is  violated,  the  QFT  compensator 
is  adjusted  and  the  time  simulation  is  repeated. 


4.1  Longitudinal  Channel  Inner  Loop  Design 

Because  the  F-16  is  statically  unstable,  the  first  step  in  the  longitudinal  channel  design  is  to 
synthesize  an  inner  feedback  loop  to  stabilize  any  plants  that  have  an  unstable  short  period  mode. 
The  feedback  parameter  used  is  pitch  rate  q.  This  step  is  typical  in  classical  flight  contol  system 
design  and  was  also  used  in  Reynold’s  design  [3,8j. 

The  overriding  priority  of  the  inner  loop  is  to  stabilize  all  of  the  plants  so  that  an  outer  loop 
can  be  accomplished  around  stable  effective  plants.  Also,  since  the  short  period  approximation  of 
the  longitudinal  dynamics  is  used,  only  the  short  period  mode  instabilities,  and  not  those  of  the 
phugoid  mode,  are  reflected. 

The  QFTCAD  is  used  to  design  this  loop  using  a  10  degree  phase  margin  angle  stability 
specification.  This  stability  specification  is  input  into  the  QFTCAD  along  with  the  structure 
shown  in  Figure  4.1.  The  previously  generated  transfer  functions  of  used  to  validate  the 
parameter  space  are  utilized  (see  Section  3.5).  The  resulting  second-order  compensator  needed  to 
stabilize  all  of  the  plants  is  given  by  Equation  (4.1). 


0.75(s-H0)(s-t-30) 

- - ^(JT^o) - 


(4.1) 
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Figure  4.1  Inner  Loop  Design  Structure 


Figure  4.2  QFTCAD  Stability  Validation  for  q  Feedback  Inner  Loop  with 

^  0.75(j  +  10)(3+30) 

—  «(5+60) 


The  corresponding  QFTCAD  stability  validation  Nichols  plot  is  depicted  in  Figure  4.2.  It  is  clear 


from  this  figure  that  all  of  the  plants  cross  the  0  dB  line  at  a  phase  angle  greater  than  -170  degrees, 


thus  satisfying  the  10  degree  phase  margin  angle  requirement. 


4.S  Longitudinal  Channel  Outer  Loop  Design 

Before  proceeding  with  an  outer  loop  design,  two  issues  must  be  addressed.  The  first  is 
selecting  an  outer  loop  feedback  parameter  and  the  second  is  validating  the  anticipated  benefits  of 
an  inner  loop  in  the  design. 
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4-8.1  Longitudinal  Outer  Loop  Control  Parameter.  Although  in  Reynold’s  thesis  the 
longitudinal  outer  loop  parameter  utilized  was  a  for  g  <  130  psf  and  C*  or  Nz  for  q  >  130  psf,  a 
different  approach  is  attempted  in  this  thesis  [8].  This  new  approach  involves  a  concerted  effort 
to  imbed  flying  qualities  into  the  design  by  selecting  an  outer  loop  parameter  which  corresponds 
closely  to  what  the  pilot  actually  desires  to  control  when  he  actuates  the  control  stick. 

Two  assumptions  are  central  in  the  development  of  the  outer  loop  parameter.  First,  that 
the  pilot  possesses  a  general  awareness  of  the  energy  state  of  the  aircraft,  and  second,  that  the 
parameter  the  pilot  wishes  to  control  varies  with  the  aircraft  energy  state.  The  author’s  expertise 
as  an  experienced  F-16  pilot  is  relied  upon  to  generate  this  dynamic  parameter. 

At  lower  energy  states,  where  little  aircraft  g  is  available,  the  pilot  commands  primarily  pitch 
rate  q  to  accomplish  his  flight  tasks.  As  the  aircraft  energy  state  increases  and  more  g  is  available, 
the  pilot  turns  more  of  his  attention  towards  managing  the  g  load,  and  less  to  managing  q.  The 
outer  loop  parameter  is  labelled  C*  and  is  referred  to  as  such  in  the  remainder  of  this  thesis.  A 
physical  explanation  for  the  proportion  of  q  and  gpu  used  to  form  C*  is  that  very  few  inflight 
tasks  demanding  precise  pitch  rate  control  by  the  pilot  are  performed  at  high  dynamic  pressures; 
however,  precise  tracking  tasks  are  common  at  low  and  medium  dynamic  pressures.  In  addition, 
at  high  energy  states  where  maximum  aircraft  g  is  available  and  sustainable,  it  is  important  that 
the  pilot  be  able  to  accurately  control  g  load  so  that  his  physical  limitations  are  not  exceeded. 
Figure  4.3  displays  the  proportion  of  q  and  gpn  that  blend  to  create  C*  at  all  subsonic  q,  where 
gpii  is  defined  as  the  positive  g  force  felt  by  the  pilot.  The  crossover  of  q  and  gpn  on  the  figure 
corresponds  approximately  to  the  dynamic  pressure  at  which  the  maximum  allowable  g  is  available 
to  the  pilot. 

The  C*  parameter  is  generated  from  the  short  period  longitudinal  state  space  plant  of  Equa¬ 
tion  (3.3)  as  follows.  First,  the  longitudinal  state  space  output  relation  y  =  Cx  -f  D«  of  Equa¬ 
tion  (3.4)  is  augmented  as  shown  in  Equation  (4.2),  where  the  shper  subscript  refers  to  the  state 
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Figure  4.3  Allocation  of  q  and  gpu  in  C*  Outer  Loop  Parameter  versus  q 


space  matrices  of  Equation  (3.3). 


y  = 


Dshper 

Bshper 


^elev 


(4.2) 


The  resulting  output  vector  of  the  augmented  longitudinal  system  is, 


a 


y  = 


9 


or 

9 


(4.3) 


The  parameter  gpu  is  now  generated  using  the  relation 


9pii  =  ~N^  =  -  9)  -  /x9] 


(4.4) 
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U  q 

(a)  No  Inner  Loc^ 


(b)  Inner  Loop  Inccaporated 

Figure  4.4  The  C*  Plant  Structures  with  and  without  an  Inner  Loop 


where  U  is  trim  velocity  in  the  x  body  axis  direction  and  is  the  distance  from  the  aircraft  center 
of  gravity  to  the  pilot  station.  Although  the  position  of  the  aircraft  center  of  gravity  is  slightly 
different  for  each  aircraft  configuration,  a  nominal  value  of  13.95  feet  is  used  to  generate  gpu.  A 
detailed  derivation  of  Equation  (4.4)  can  be  found  in  Reynold’s  thesis  [8].  Finally,  q  and  Qpu  are 
blended  to  produce  the  dynamic  pressure  dependent  C*  parameter  profile  depicted  in  Figure  4.3. 

4-2.2  Inner  Loop  Benefits.  To  evaluate  the  benefit  of  an  inner  loop  in  this  design, 
frequency  templates  are  generated  for  C*  with  and  without  an  inner  loop.  If  the  QFT  frequency 
templates  for  the  system  with  an  inner  loop  are  smaller  in  the  desired  frequency  range,  the  inner 
loop  is  deemed  to  be  of  benefit;  otherwise,  the  incorporation  of  an  inner  loop  not  only  complicates 
the  design  but  results  in  an  outer  loop  compensator  of  higher  order. 

The  necessary  plant  transfer  functions  for  C*  are  generated  by  Matlab  using  the  structures 
illustrated  in  Figure  4.4  and  then  loaded  into  the  QFTCAD. 
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Figure  4.5  QFTCAD  Frequency  Templates  for  System  with  Inner  Loop 

The  QFTCAD  frequency  templates  for  the  longitudinal  system  with  an  inner  loop  are  shown 
in  Figure  4.5.  These  templates  are  compared  to  those  in  Figure  4.6,  which  are  for  the  system 
without  an  inner  loop.  By  examining  these  two  figures,  it  is  noted  that  the  inner  loop  feedback 
does  dramatically  decrease  the  template  size  at  low  frequencies,  but  it  increases  the  template 
size  in  the  middle  and  high  frequencies  of  interest  in  the  design.  Uncertainty  at  low  frequencies 
is  much  easier  to  compensate  for  than  uncertainty  in  the  middle  to  high  design  frequencies.  In 
addition,  uncertainty  below  the  typical  control  bandwidth  of  a  pilot  can  be  essentially  dismissed. 
Consequently,  the  inner  loop  does  not  ease,  but  hinders  the  design  effort  and  is  therefore  not 
incorporated. 

4-2.3  Longitudinal  QFT  Specifications.  With  the  outer  loop  parameter  and  inner  loop 
benefit  issues  settled,  the  next  and  necessary  step  in  the  design  of  the  outer  loop  is  the  definition 
of  the  specification  models  required  as  input  to  the  QFTCAD.  The  document  used  as  the  primary 
source  for  flying  quality  specifications  is  Military  Standard  1797A  (Mil  Std)  [2].  The  specifications 
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Figure  4.6  QFTCAD  Frequency  Templates  for  System  without  Inner  Loop 


IIK^ 

At 

Table  4.1  Mil  Std  Longitudinal  Time  Domain  Tracking  Specifications 

from  the  Mil  Std  used  in  this  design  are  those  necessary  for  the  achievement  of  Level  1  flying 
qualities. 


4-2.3. 1  Tracking  Specifications.  In  the  Mil  Std  the  time  domain  tracking  specifica¬ 
tions  are  defined  by  analyzing  a  typical  pitch  rate  response  to  a  stick  force  step  input.  Figure  4.7 
depicts  a  typical  underdamped  response  to  a  step  input  with  the  necessary  parameters  used  to  de¬ 
fine  time  domain  tracking  specifications  labelled.  Table  4.1  sets  forth  the  specifications  needed  to 
achieve  Level  1  flying  qualities,  where  V  is  the  aircraft  velocity  in  feet  per  second.  The  ti  parameter 
defines  a  maximum  onset  delay  requirement,  the  At  parameter  dictates  an  onset  rate  requirement, 
and  the  ratio  imposes  a  minimum  damping  requirement. 
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Figure  4.7  Step  Response  used  to  Define  Longitudinal  Response  Specifications 


A  couple  of  difficulties  exist  when  attempting  to  map  these  time  domain  specifications  into 
the  frequency  domain.  First,  the  specifications  in  Table  4.1  are  based  only  on  pitch  rate  response, 
but  the  outer  loop  parameter  of  this  design  is  not  solely  pitch  rate.  Second,  the  allowable  range  of 
At  varies  with  flight  condition  (aircraft  velocity)  and  is  therefore  a  floating  specification. 

The  first  of  these  difficulties  is  addressed  in  the  following  way.  At  the  lower  aircraft  energy 
levels  (i.e.,  at  low  airspeeds),  C*  is  primarily  composed  of  pitch  rate.  The  onset  delay  (ti)  and  onset 
rate  {At)  requirements  are  most  difficult  to  satisfy  at  low  dynamic  pressures  due  to  reduced  control 
surface  effectiveness.  Also,  the  Mil  Std  does  not  set  forth  any  explicit  time  domain  specifications 
for  Qpii  response  to  pilot  pitch  stick  commands.  Hence,  the  pitch  rate  step  response  specifications 
are  used  as  the  starting  point  for  synthesizing  QFT  upper  and  lower  tracking  bound  models. 

Next,  the  issue  of  the  floating  At  specifications  is  settled  by  selecting  reasonable  upper  and 
lower  limits  on  At  based  on  the  fact  the  lower  velocity  plants  respond  more  slowly  to  pitch  stick 
commands  than  the  higher  velocity  plants. 
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The  initial  form  of  the  upper  tracking  bound  transfer  function  model  is 


Tru  = 


+  2Cai„  +a)„2' 


(4.5) 


It  is  found  that  ^  must  be  greater  than  or  equal  to  0.35  in  order  to  satisfy  the  damping 
specification  from  Table  4.1.  A  C  of  0-4  is  selected  for  the  upper  tracking  bound  model. 

Next,  an  cj„  of  0.7  radians  per  second  is  selected  for  two  reasons.  First,  most  of  the  plant 
uncertainty  occurs  at  frequencies  less  than  or  equal  to  1  radian  per  second  (see  templates  in  Fig¬ 
ure  4.6).  Second,  placing  the  responses  with  the  greatest  overshoot  and  slowest  settling  in  the  lower 
frequencies  allows  the  pilot  the  most  time  to  compensate  for  these  undesired  response  characteris¬ 
tics,  or  possibly  places  these  characteristics  in  frequencies  below  the  bandwidth  of  the  pilot. 

The  upper  limit  of  a  typical  pilot’s  bandwidth  is  approximately  3.5  rps  [7].  Consequently,  the 
upper  tracking  bound  model  must  have  a  bandwidth  of  at  least  3.5  radians  per  second.  A  target 
bandwidth  of  5  radians  per  second  is  chosen  for  the  upper  bound  and  is  defined  as  the  frequency 
at  which  the  magnitude  of  the  frequency  response  of  the  upper  bound  model  descends  through  -3 
dB.  It  is  necessary  to  add  a  zero  to  the  model  in  order  to  achieve  the  required  bandwidth. 

The  time  response  overshoot  Mo  =  Agi  is  not  given  a  maximum  value  by  the  Mil  Std.  A  value 
of  Mo  =  0.75  is  arbitrarily  chosen  as  this  maximum  value  to  allow  for  a  small  At.  A  pole-zero  pair 
is  added  to  the  model  to  achieve  these  characteristics  while  complying  with  the  target  bandwidth 
previously  set  forth.  The  resulting  upper  tracking  bound  model  is 


Tru 


5(s  +  0.7)2 

(s2-f0.56s-f  0.49)(s  +  5)’ 


(4.6) 
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Model 

mm 

At 

Tru 

0 

0.215 

0.255 

Trl 

0.113 

1.132 

overdamped 

Table  4.2  QFT  Upper  and  Lower  Tracking  Bound  Model  Step  Response  Characteristics 

The  primary  concern  in  the  synthesis  of  the  lower  tracking  bound  model  is  that  the  onset 
rate  ti  does  not  exceed  0.12  seconds.  The  initial  model  is 


Trl  = 


2 

s  "L  2 


(4.7) 


chosen  to  ensure  a  bandwidth  of  2  radians  per  second.  The  model  is  then  augmented  with  another 
pole  to  fulfill  the  QFT  requirement  that  the  upper  and  lower  bounds  diverge  above  the  frequency 
the  upper  bound  crosses  through  0  dB.  This  additional  pole  is  placed  at  the  lowest  frequency 
possible  without  increasing  the  onset  rate  beyond  0.12  seconds.  The  resulting  lower  tracking  bound 
model  is 

^“=(.  +  2K»  +  3)'  (“> 

These  upper  and  lower  QFT  tracking  bound  models  require  all  plants  to  have  an  onset  rate 
less  than  0.12  seconds,  M,  <  0.75,  a  ratio  less  than  0.3,  and  a  bandwidth  between  2  and  5 
radians  per  second.  Table  4.2  lists  the  exact  characteristics  of  the  step  input  time  responses  of  the 
tracking  models.  In  addition.  Figure  4.8  depicts  both  the  frequency  response  magnitude  and  step 
input  time  response  plots  of  the  QFT  tracking  bound  models. 

4-2.3. 2  QFT  Siability  Specifications.  The  required  stability  specifications  are  clearly 
set  forth  in  the  Mil  Std.  The  open- loop  frequency  response  of  all  plants  must  have  a  phase  margin 
angle  of  at  least  30  degrees  and  a  gain  margin  of  at  least  6  dB.  The  30  degree  phase  margin  angle 
requirement  is  entered  directly  into  the  QFTCAD  to  generate  the  stability  bounds,  and  the  6  dB 
gain  margin  is  verified  visually  on  the  QFTCAD  Nichols  chart  displays. 
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Figure  4.8  Frequency  Response  Magnitude  (top)  and  Step  Input  Time  Response  (bottom)  of 
QFT  Upper  and  Lower  Longitudinal  Tracking  Bound  Models 

4.2.S.3  QFT  Disturbance  Specifications.  A  thorough  discussion  of  modelling  eleva¬ 
tor  trim  changes  as  an  external  disturbance  is  presented  in  Reynold’s  thesis  [8].  It  is  uncertain, 
however,  as  to  what  level  of  disturbance  rejection  is  desired  or  required  in  order  to  prevent  the 
elevator  trim  changes  from  posing  a  problem  to  the  pilot  during  longitudinal  tracking  tasks.  In 
light  of  this  uncertainty  and  considering  the  difficulty  of  satisfying  even  the  tracking  specifications 
required  in  this  design,  disturbance  rejection  is  not  prioritized  above  tracking.  Instead,  disturbance 
rejection  is  treated  as  a  by  product  of  meeting  tracking  specifications  as  opposed  to  being  pursued 
as  an  end  in  itself.  This  is  not  unreasonable  considering  the  fact  that  the  longitudinal  compensator 
of  this  design  is  required  to  have  pure  integral  action  which  gives  the  system  inherent  disturbance 
rejection  properties. 

4.2.4  F-16  Longitudinal  Performance  Requirements.  To  achieve  performance  compara¬ 

ble  to  the  current  F-16  FCS,  performance  benchmarks  must  be  established.  The  performance 
benchmarks  used  in  this  design  are 

1.  a  -  attain  25®,  but  do  not  exceed  30® 
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2.  Qpn  -  attain  9  g's,  but  do  not  exceed  9.45  g’s 


Only  the  benchmark  which  is  achieved  first  at  each  dynamic  pressure  is  demanded  of  the  FCS.  For 
instance,  at  low  dynamic  pressures  only  benchmark  1  applies  and  at  high  dynamic  pressures  only 
benchmark  2  applies. 

4-2.5  QFT  Longitudinal  Channel  Design.  The  plan  of  attack  for  the  outer  loop  design  is 
to  synthesize  the  simplest  compensator  possible  (maximum  of  third-order  and  with  poles  <  60  rps) 
to  meet  the  tracking  specifications  for  the  frequency  range  of  0.5  to  3.5  radians  per  second.  The 
design  is  also  required  to  meet  stability  specifications  for  all  frequencies.  If  a  single  compensator 
cannot  achieve  the  tracking  or  stability  requirements,  then  either  the  gain  of  the  compensator  will 
be  scheduled  on  q,  or  a  completely  different  compensator  will  be  generated  to  be  switched  to  at  a 
particular  q. 

4-2.5. 1  QFTCAD  Plant  Input.  The  C*  plants  previously  generated  by  Matlab  from 
the  structure  of  Figure  4.4(a)  are  used  in  the  outer  loop  design.  Once  the  plants  are  loaded,  the 
effective  plants  P*  (which  include  the  actuator  dynamics)  and  the  Q  matrices  are  generated.  In 
this  generation  process,  the  QFTCAD  cancels  poles  and  zeros  within  0.001  of  one  another  and 
eliminates  all  poles  and  zeros  outside  of  the  dynamic  frequency  range  of  10“®  to  10“*  rps.  The 
resulting  plant,  effective  plant,  and  q  transfer  functions  for  all  plants  are  listed  in  the  supplement 
to  this  thesis. 

4-2. 5.2  QFTCAD  Structure.  The  structure  entered  into  the  QFTCAD  for  the 
longitudinal  design  is  portrayed  in  Figure  4.9.  Recall  that  the  proportion  of  pitch  rate  and  g^n 
which  make  up  C*  is  dependent  on  the  current  dynamic  pressure  of  the  aircraft  (see  Figure  4.3). 
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Figure  4.9  QFTCAD  Longitudinal  Channel  Structure 


4. 2. 5. 3  QFT  Design  Boundaries.  The  specifications  discussed  in  Section  4.2.3  are 
entered  into  the  QFTCAD  to  be  used  to  generate  design  boundaries  on  the  Nichols  chart.  However, 
before  the  boundaries  can  be  generated,  two  items  must  be  accomplished. 

First,  frequency  templates  must  be  formed.  Recall  that  in  this  design  templates  are  formed 
for  the  frequency  spectrum 


w  =  {0.05, 0.075, 0.1 , 0.25, 0.5, 0.75, 1, 3, 5, 7, 10, 15,20, 25, 30, 50}  rps. 


(4.9) 


Second,  a  nominal  plant  is  selected.  Plant  281  is  selected  as  the  nominal  plant  due  to  its 
position  at  the  top  of  the  30  rps  template.  The  frequency  templates  for  the  outer  loop  are  shown 
in  Figure  4.6  with  plant  281  shown  as  the  nominal. 

The  QFTCAD  now  generates  the  tracking  and  stability  boundaries  to  be  used  in  the  design. 
External  disturbance  boundaries  are  not  generated  because  of  the  reasons  discussed  in  Sections  3.2 
and  4. 2. 3. 3. 

4-2.5. 4  Ouier  Loop  Compensator  Design.  The  QFTCAD  ‘Design  Compensator’ 
feature  is  selected,  and  the  nominal  open-loop  function  is  shaped  through  the  synthesis  of  the 
outer  loop  compensator  Gc*  •  To  ensure  zero  steady-state  step  input  tracking  error  and  to  enhance 
disturbance  rejection,  Gc-  is  given  a  pole  at  the  origin  to  make  the  forward  transfer  function 
Type  1  [5].  Next,  poles  and  zeros  are  added  to  the  compensator  so  that  the  nominal  loop  satisfies 
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Figure  4.10  QFTCAD  Longitudinal  Frequency  Templates  for  Clean,  One  Fuel  Tank,  and  Two 
Fuel  Tank  Configurations 


the  necessary  tracking  and  stability  specifications  set  forth  in  Section  4.2.3.  The  compensator 
is  then  simulated  using  Mailab  to  observe  the  time  domain  performance  of  all  of  the  plants  and 
to  identify  any  actuator  rate  or  control  surface  deflection  saturations.  During  this  process,  the 
first-order  prefilter 


Fc- 


3.5 

5  -1-  3.5 


(4.10) 


is  used  due  to  its  bandwidth  of  3.5  rps. 


In  the  first  design  attempt,  a  single  compensator  of  third-order  could  not  shape  the  nominal 
loop  to  satisfy  all  of  the  tracking  and  stability  boundaries.  Rather  than  designing  two  separate 
compensators  or  applying  gain  scheduling,  it  is  decided  to  eliminate  all  of  the  three  fuel  tank 
configured  plants  and  then  reaccomplish  the  design.  The  three  tank  F-16  configuration  is  normally 
used  only  in  long  range  deployment  missions,  and  not  in  weapon  employment  missions.  Therefore, 
Level  1  flying  qualities  are  required  of  only  the  clean,  one  fuel  tank,  and  two  fuel  tank  configurations. 
The  QFTCAD  frequency  templates  for  the  remaining  199  plants  are  displayed  in  Figure  4.10. 
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Figure  4.11  QFTCAD  Nominal  Loop  Shaping  Nichols  Chart  for  Gc-  ~  ^ 


After  removing  the  three  tank  configurations,  the  QFTCAD/time  simulation  process  de¬ 
scribed  in  the  previous  paragraphs  results  in  the  compensator 


Gc» 


3.3(s+  1.7)(s -1-3.5) 
s(s  +  12) 


(4.11) 


The  QFTCAD  nominal  loop  shape  along  with  tracking  and  stability  bounds  is  depicted  on  the 
Nichols  chart  in  Figure  4.11.  This  figure  clearly  shows  that  the  nominal  loop  has  an  of  approxi¬ 
mately  20  rps  and  a  gain  margin  of  greater  than  6  dB.  Since  the  nominal  plant  is  located  at  the  top 
of  the  frequency  templates  from  20  to  50  rps,  all  of  the  other  plant  loops  have  an  less  than  and 
a  gain  margin  greater  than  the  nominal  loop.  This  figure  also  shows  that  the  very  low  frequency 
tracking  bounds  are  violated,  but  the  bounds  in  the  desired  range  of  0.5  to  3.5  rps  are  satisfied. 

The  ‘Design  Prefilter’  QFTCAD  display  for  Fc»  =  is  shown  in  Figure  4.12.  Since  the 
very  low  frequency  tracking  bounds  are  not  satisfied  during  the  loop  shaping,  the  upper  and  lower 
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Figure  4.12  QFTCAD  Prefilter  Nichols  Chart  for  Fc>  = 

prefilter  bounds  cannot  be  satisfied  at  these  frequencies,  but  in  the  frequency  range  of  0.5  to  10 
rps,  the  prefilter  bounds  are  satisfied. 

4.2.6  Longitudinal  Design  Time  Domain  Analysis.  The  first  step  in  the  time  domain 
analysis  is  to  subject  the  compensated  system  to  a  unit  step  input.  The  results  of  this 

simulation  for  a  representative  sample  of  plants  spanning  the  full  range  of  q  are  illustrated  in 
Figures  4.13  and  4.14.  The  plots  of  Seiev  and  a  depict  perturbations  from  the  initial  condition  trim 
values,  and  not  the  total  values  including  trim.  The  ballooning  above  1 .25  observed  in  the  C*  step 
responses  is  for  aircraft  plants  corresponding  to  dynamic  pressures  below  200  pounds  per  square 
foot  (psf).  The  q  response  of  every  plant  is  evaluated  against  the  requirements  of  Table  4.1.  A 
tabulation  of  the  results  of  this  evaluation  is  contained  in  Appendix  C.l.  An  inspection  of  the  data 
in  the  appendix  confirms  that  the  onset  delay,  onset  rate,  and  damping  requirements  of  Table  4.1 
are  met  by  all  plant  models. 

The  method  applied  in  this  design  to  avoid  physical  saturations  while  meeting  performance 
requirements  is  to  scale  the  stick  force  input  of  the  pilot  rather  than  to  directly  schedule  the 
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Figure  4.14  Compensated  System  Time  Response  to  Unit  Step  Input  (2  of  2) 


compensator  gain.  To  accomplish  the  stick  force  scaling,  a  maximum  profile  scheduled  on  q  is 
synthesized  using  the  results  of  the  unit  step  input  time  simulation.  The  data  generated  for 
each  plant  by  Equations  (4.12)  through  (4.15)  is  plotted  versus  q  to  form  upper  and  lower  bounds 
for  this  profile. 
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The  maximum  at  a  particular  q  must  lie  below  the  plot  of  Equations  (4.12)  and  (4.13)  in  order 
to  avoid  elevator  deflection  and  actuator  rate  limitations,  respectively.  Also,  to  achieve  performance 
comparable  to  the  current  F-16  FCS,  the  must,  for  all  q,  allow  the  aircraft  to  reach  either 

25  degrees  of  a  or  9  g^s  at  the  pilot  station,  whichever  occurs  first.  This  condition  is  met  when  the 
C'rmd  profile  is  above  the  lowest  of  the  two  plots  of  Equations  (4.14)  and  (4.15).  The  resulting 
maximum  profile  is  displayed  in  Figure  4.15  along  with  the  upper  and  lower  bounds.  In  order 
to  help  declutter  the  presentation,  the  lower  bound  reflecting  desired  a  and  gpu  performance  is 
combined  into  a  single  bound  by  plotting  only  the  minimum  of  the  two  individual  bounds. 


Next,  the  longitudinal  FCS  is  subjected  to  step  inputs  equal  to  the  maximum  from 
Figure  4.15  to  validate  the  expected  performance.  The  results  of  this  simulation  are  recorded  in 
Figures  4.16  and  4.17.  It  is  evident  in  these  figures  that  all  of  the  aircraft  plants  not  only  meet, 
but  some  exceed  the  desired  a  and  gpn  performance.  A  listing  of  the  results  of  this  simulation  for 
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Figure  4.15  Maximum  Profile  vs.  q  Along  with  Upper  and  Lower  Bounds 

all  plants  is  contained  in  Appendix  C.2.  The  fact  that  the  desired  maximum  values  of  a  and  gpn 
are  exceeded  warrants  the  discussion  of  limiters  in  the  next  section. 

4. 2. 6.1  Longitudinal  Channel  Limiters.  Due  to  concern  about  loss  of  aircraft  con¬ 
trol,  loss  of  pilot  consciousness,  or  overstress  of  the  aircraft  itself,  the  need  for  longitudinal  channel 
limiters  is  now  assessed.  In  order  to  ensure  aircraft  stability,  a  maximum  transient  a  of  30  degrees 
is  desired.  Also,  a  meiximum  overshoot  of  five  percent  (0.45  g  for  a.9  g  aircraft)  in  gpu  is  pursued. 
Finally,  an  upper  limit  of  6  g’s  per  second  is  placed  on  g  onset  rate  in  order  to  reduce  the  likeli¬ 
hood  of  pilot  loss-of-consciousness  during  the  initiation  of  high  g  maneuvers.  The  time  responses 
of  Figures  4.16  and  4.17  indicate  that  limiters  are  required  to  satisfy  these  requirements. 

An  ad  hoc  limiter  scheme  involving  rate  and  magnitude  limiting  of  a  and  gpu  is  concocted 
to  satisfy  the  requirements  set  forth  in  the  previous  paragraph.  Since  a  and  gpu  are  not  each 
independently  controlled  by  the  limiting  scheme  is  not  as  straight  forward  as  simple  regulation 

through  feedback.  Instead,  the  input  by  the  pilot  is  reduced  through  feedback  as  rates  or 

magnitudes  exceed  certain  values.  Figure  4.18  graphically  illustrates  the  ad  hoc  limiters  used  in 
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Figure  4.17  Compensated  System  Time  Response  to  Maximum  Step  Input  (2  of  2) 
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this  design.  The  graphs  indicate  how  is  scaled  according  to  the  rates  and  magnitudes  of  a 

and  gpii-  Plots  (a)  and  (c)  indicate  that  no  magnitude  limiting  is  initiated  until  above  the  desired 
ffpii  and  a  performance  requirements  of  9  g’s  and  25  degrees,  respectively.  Any  limiting  below  these 
values  is  the  result  of  rate  limiting.  The  action  of  the  g  rate  limiter  in  Figure  4.18(b)  during  the 
initiation  of  a  high  g  maneuver  is  to  keep  the  g  onset  rate  from  exceeding  the  6  g  per  second  limit 
early  in  the  maneuver,  to  stabilize  this  onset  rate,  and  then  to  slow  the  onset  rate  again  when 
approaching  the  maximum  g  of  9.45  to  prevent  overshoot.  The  a  rate  limiter  of  plot  (d)  slows  the 
a  rate  so  as  to  prevent  an  overshoot  of  the  30  degree  maximum.  The  two  dimensional  look-up 
tables  used  in  Matlab  to  simulate  the  rate  limiters  of  plots  (b)  and  (d)  are  found  in  Appendix  D. 
Figure  4.22  indicates  how  the  limiters  fit  into  the  design  structure. 

4- £.7  Longitudinal  Design  Validation.  The  QFTCAD  is  used  to  perform  frequency  domain 
validation  and  Matlab  is  used  to  perform  time  domain  validation  of  the  design. 

4-2.7. 1  QFTCAD  Frequency  Domain  Validation.  To  confirm  that  all  of  the  remain¬ 
ing  198  open-loop  transfer  functions  meet  the  30  degree  phase  margin  angle  and  6  dB  gain  margin 
stability  requirements,  the  QFTCAD  ‘Stability  Validation’  feature  is  selected.  Figure  4.19  depicts 
the  resulting  Nichols  chart.  Examination  of  this  figure  reveals  that  all  199  open-loop  functions  do 
in  fact  meet  the  desired  stability  specifications. 

Next,  by  selecting  the  QFTCAD  ‘Tracking  Validation’  feature,  the  closed-loop  frequency 
responses  for  all  199  plants  are  plotted  along  with  the  upper  and  lower  frequency  domain  tracking 
bounds,  Tru  and  Trl.  These  Bode  plots  are  displayed  in  Figure  4.20.  Both  the  anticipated 
violation  of  the  tracking  bounds  at  the  very  low  frequencies  and  the  satisfaction  of  the  tracking 
bounds  between  0.5  and  3.5  rps  are  evident;  hence,  the  design  achieves  the  desired  frequency  domain 
tracking  specifications. 
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Percentage  C*(cmd)  Percentage  C*(cmd) 


(a)  g  Limiting  (b)  g  Onset  Rate  Limiting 


alpha  (deg) 
(c)  alpha  Limiting 


alpha  (deg)  alpha  rate  (deg/sec) 

(d)  alpha  Rate  Limiting 


Figure  4.18  Graphical  Display  of  a  and  gpn  Limiters 
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Figure  4.21  QFTCAD  Tracking  Validation  Bode  Plots 


The  QFTCAD  ‘External  Disturbance  Rejection  Validation ’feature  generates  Bode  magnitude 
plots  for  each  plant  indicating  the  level  of  external  disturbance  rejection  achieved.  The  resulting 
Bode  plots  for  the  longitudinal  system  are  shown  in  Figure  4.21  where  the  external  disturbance 
considered  is  the  change  in  elevator  trim  as  a  result  of  changes  in  flight  condition.  Reynolds 
estimated  that  for  w  <  0.157  rps  the  effects  of  elevator  trim  disturbances  should  be  considered  [8]. 
Applying  this  frequency  range  estimate  to  the  plots  of  Figure  4.21,  it  is  concluded  that  a  respectable 
disturbance  rejection  of  -20  dB  is  achieved  by  the  longitudinal  design. 

4.S.7.2  Longitudinal  Design  Time  Domain  Validation.  A  representative  set  of  sys¬ 
tem  time  responses  to  maximum  step  inputs  with  the  longitudinal  limiter  scheme  in  place  is 
displayed  in  Figures  4.23  and  4.24.  A  comparison  of  these  figures  to  Figures  4.16  and  4.17  shows 
that  the  limiter  scheme  accomplishes  its  intended  purpose.  A  tabular  listing  of  the  results  of  this 
simulation  for  all  plants  with  the  limiters  in  place  is  found  in  Appendix  C.3.  These  figures  and 
the  appendix  data  demonstrate  that  the  design  meets  the  performance  benchmarks  while  avoiding 
actuator  rate  and  control  surface  deflection  saturations. 
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Figure  4.22  Longitudinal  Limiter  Structure 


4.3  Final  Longitudinal  Channel  Design 

The  final  design  in  the  longitudinal  channel  is  best  described  by  the  block  diagram  illustrated 
in  Figure  4.25.  To  summari2e,  the  pitch  stick  force  from  the  pilot  is  first  scaled  by  the  maximum 
C*^d  profile  and  then  scaled  by  the  limiters.  This  scaled  pitch  stick  force  is  then  filtered  and 
enters  the  C*  compensator  feedback  loop,  where  it  is  transformed  into  an  appropriate  input  to  the 
elevator.  Both  the  prefilter  and  compensator  of  the  longitudinal  design  satisfy  the  requirements  set 
forth  in  Section  2.3. 
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Figure  4.23 


Compensated  System  Time  Response  to  Maximum  Ccm.  Step  Input  with  Limiters 
in  Place  (1  of  2) 
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Figure  4.25  Final  Longitudinal  FCS  Design  Block  Diagram 
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V.  Lateral  Channel  Design 


The  lateral  design  follows  the  same  integrated  approach  that  is  used  in  the  longitudinal 
channel  design  (see  introduction  to  Chapter  IV). 


5.1  Dutch  Roll  Damping  Loop  Design 

In  line  with  traditional  lateral  FCS  design  and  Reynold’s  research,  the  first  step  taken  is  to 
improve  the  dutch  roll  mode  damping  of  the  basic  lateral  aircraft  plants  [3,8].  This  is  done  by 
feeding  yaw  rate  r  back  through  a  washout  filter  of  the  form 


to  the  rudder.  This  loop  acts  to  limit  undesired  yaw  rates  caused  by  flaperon  and  horizontal  tail 
deflection  during  the  initiation  of  rolling  maneuvers  for  a  transient  period  equal  to  r  seconds,  while 
not  interfering  with  the  steady  state  yaw  rates  of  a  sustained  turning  maneuver  [3].  A  block  diagram 
of  this  dutch  roll  damping  loop  is  illustrated  in  Figure  5.1. 

A  root  locus  analysis  of  the  open-loop  transfer  functions  (OLTF’s)  from  rudder  input  to  yaw 
rate  output  is  performed.  The  transfer  functions  of  the  lateral  aircraft  plants  corresponding  to  the 
maximum  and  minimum  dynamic  pressure  are  chosen  for  this  analysis.  The  objective  of  the  analysis 
is  to  determine  values  of  K  and  r  of  the  washout  filter  described  by  Equation  (5.1)  that  achieve  the 

^rud(cmd) 
reference 
signal 


Figure  5.1  Dutch  Roll  Damping  Washout  Circuit  Block  Diagram 


5-1 


(ft)  Low  Dynamic  Pretsurt  Plant 


Real  Axis 


Figure  5.2  Dutch  Roll  Damping  Loop  Root  Locus  Plots  for  t  —  3.33  seconds 

best  damping  ratio  (^)  of  the  dutch  roll  mode  poles  of  both  the  high  and  low  q  closed-loop  transfer 
functions  (CLTF’s). 

Although  the  C  of  the  high  q  plant  is  improved  dramatically  by  the  introduction  of  the  washout 
filter,  the  damping  ratio  of  the  low  q  plant  is  not.  Consequently,  the  r  selected  is  the  value  which 
results  in  the  greatest  possible  damping  ratio  for  the  low  q  plant  root  locus.  The  r  that  accomplishes 
this  is  3.33  seconds,  and  is  the  same  eis  that  used  by  Reynolds  [8].  The  resulting  root  locus  plots 
for  the  high  and  low  q  plant  OLTF’s  are  depicted  in  Figure  5.2. 

Next,  a  value  of  K  is  selected  to  maximize  the  damping  of  the  closed-loop  dutch  roll  mode 
poles  of  both  of  the  loci  of  Figure  5.2.  Although  a  if  of  4  maximizes  the  C  of  the  low  q  plant  dutch 
roll  mode  poles,  it  pushes  the  high  q  plant  dutch  roll  mode  towards  instability.  Hence,  a  smaller 
gain  must  be  chosen  as  a  compromise.  A  gain  of  2.6  achieves  essentially  equal  closed-loop  dutch 
roll  mode  damping  for  both  plants.  The  -f’s  in  Figure  5.3  show  the  position  of  the  closed-loop 
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(•)  Low  Dyiwne  Preuurt  Rinl 


(b)  High  Dyntn^c  Pretiure  P1«l 


Figure  5.3  Root  Locus  Plots  and  Closed  Loop  Roots  for  Hu>o  = 
dutch  roll  mode  roots  of  the  high  and  low  q  plants  for  the  washout  filter  of  Equation  (5.2). 


HjifO  — 


2.6s 
s  +  0.3 


(5.2) 


The  right  is  reserved  to  adjust  K  once  the  outer  loops  are  in  place,  since  the  outer  loops  will  no 
doubt  interact  to  some  extent  with  the  dutch  roll  mode  damping  loop. 


5.S  Lateral  Channel  Outer  Loop  Design 

5.2.1  Lateral  Outer  Loop  Control  Parameters.  As  in  Reynold’s  thesis,  the  lateral  outer 
loop  control  parameters  in  this  design  are  roll  rate  p  and  sideslip  angle  0  [8].  In  this  design,  the 
emphasis  is  placed  on  the  p  feedback  control  loop  since  p  commands  are  rarely  used  by  pilots  to 
accomplish  inflight  lateral  tracking  tasks. 
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5.S.2  Lateral  QFT  Specifications.  As  in  the  longitudinal  channel,  the  Military  Standard 
1797A  is  used  as  the  basis  for  the  synthesis  of  frequency  domain  specifications  to  use  in  the  QFT 
portion  of  the  lateral  design  [2]. 


5. 2. 2.1  QFT  Tracking  Specifications.  In  the  roll  channel,  the  Mil  Std  specifies  a 
minimum  roll  mode  time  constant  of  Troit  =  1  second.  This  is  interpreted  as  a  4  second  maximum 
limit  on  step  response  settling  time.  The  Mil  Std  also  requires  that  the  aircraft  be  able  to  roll 
through  360  degrees  of  bank  in  2.8  seconds  or  less  in  order  to  achieve  Level  1  flying  qualities. 
An  additional  Level  1  combat  specification  is  the  capability  to  roll  through  90  degrees  of  bank 
in  1  second  or  less.  These  bank  angle  change  specifications  are  interpreted  as  applicable  only  to 
reasonable  aircraft  energy  levels,  and  not  to  flight  conditions  of  low  dynamic  pressure. 

The  Mil  Std  also  stipulates  limits  on  the  cross-coupling  of  roll  rate  commands  into  the  yaw 
channel.  It  states  that  a  roll  rate  command  of  1  degree  per  second  result  in  no  more  than  0.022 
degrees  of  sideslip  except  at  low  airspeeds  where  it  must  be  held  to  a  maximum  of  0.067  degrees. 
Also,  the  sideslip  resulting  from  a  maximum  roll  rate  command  must  not  exceed  6  degrees. 


Roll  Rate  QFT  Tracking  Specification  Models.  The  only  one  of  the  Mil  Std 
specifications  used  in  the  synthesis  of  QFT  upper  and  lower  roll  rate  response  tracking  models  is 
the  roll  mode  time  constant  minimum  of  1  second.  Although  time  to  roll  performance  requirements 
are  set  forth  by  the  Mil  Std,  they  are  not  used  in  the  development  of  tracking  specification  models. 
Instead,  these  requirements  are  examined  during  time  domain  simulation. 


The  upper  p  tracking  specification  model  is  chosen  as  an  underdamped  response  with  a  ^ 
of  0.5  whose  settling  time  does  not  exceed  4  seconds.  By  applying  the  two  percent  settling  time 
formula 


T,  = 


4 


CWn’ 


(5.3) 
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Figure  5.4 


Roll  Rate  QFT  Upper  and  Lower  Tracking  Specification  Model  Frequency  and  Step 
Time  Responses 


a  settling  time  of  4  seconds  dictates  an  undamped  natural  frequency  of  w„  =  2  rps  [5].  The  Tru 
model  includes  a  zero  to  achieve  divergence  with  Tri  above  Tru's  0  dB  crossing  frequency.  This 
zero  is  also  placed  to  set  the  model’s  bandwidth  to  5  rps.  The  resulting  model  is 


Tru  = 


4(s  +  1) 

+  2s  +  4 


(5.4) 


The  lower  p  tracking  specification  model  is  initially  chosen  as  to  set  a  bandwidth  of  2  rps. 
Then,  this  initial  model  is  augmented  with  a  pole  to  promote  divergence  with  Tru  as  discussed  in 
the  last  paragraph.  This  additional  pole  is  placed  such  that  the  settling  time  of  the  model  does  not 
exceed  4  seconds.  The  resulting  model  is 


(s+1.25)(s  +  2)’ 

Figure  5.4  depicts  both  the  frequency  and  step  input  time  responses  of  Tru  and  Trr  for  the  roll 
rate  loop. 
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As  mentioned  before,  /?  is  not  typically  commanded  by  the  pilot;  hence,  the  cross-coupling 
effect  of  0ctnd  to  p  is  of  little  concern.  In  the  rare  case  that  the  pilot  does  command  sideslip, 
the  rudder  inputs  from  the  pilot  are  filtered  as  necessary  to  reduce  any  undesired  cross-coupling 
effects.  The  required  filtering  of  sideslip  commands  is  determined  from  the  results  of  time  domain 
simulation. 

On  the  other  hand,  the  cross-coupling  effects  of  rudder  deflection  as  a  result  of  the  dutch 
roll  damping  and  /3  regulation  loops  definitely  warrant  attention.  It  is  not  evident  what  level  of 
rejection  is  required  to  eliminate  any  undesirable  cross-coupling  in  the  roll  rate  response,  so  the 
roll  rate  loop  compensator  is  designed  assuming  no  cross-coupling.  The  roll  rate  command  time 
response  of  the  system  is  analyzed  to  determine  if  enough  cross-coupling  rejection  is  present.  If  the 
yaw  damping  and  /?  loops  adversely  affect  the  roll  rate  response,  the  roll  rate  compensator  and/or 
these  two  loops  resposible  for  the  cross-coupling  are  adjusted. 

Sideslip  QFT  Tracking  Specification  Models.  The  limits  imposed  by  the  Mil 
Std  on  sideslip  angle  due  to  rolling  maneuvers  is  modelled  as  an  upper  cross-coupling  bound  for 
the  S  feedback  loop  design.  A  constant  cross-coupling  upper  bound  of  0,067  is  used.  The  Mil  Std 
stipulates  no  explicit  requirements  for  the  sideslip  angle  response  to  pilot  rudder  inputs;  therefore, 
no  tracking  models  are  required  for  the  fi  response.  However,  the  QFTCAD  requires  the  input 
of  tracking  bounds  for  this  MISO  loop  in  order  to  generate  the  cross-coupling  bounds  and  to  use 
the  ‘Design  Prefilter’  feature,  so  arbitrary  tracking  models  are  selected  and  entered  into  the  CAD. 
The  p  loop  upper  tracking  bound  model  of  Equation  5.4  is  adopted  for  the  fi  loop  and  the  lower 
tracking  bound  model 

is  utilized. 
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Figure  5.5  Lateral  Channel  QFT  Tracking  Boundaries 


Lateral  Tracking  Bound  Summary.  The  tracking  bounds  for  all  four  QFT 
MISO  loops  are  illustrated  in  Figure  5.5 

5. 2. 2.2  QFT  Lateral  Stability  Specifications.  As  in  the  longitudinal  channel,  the  p 
and  /?  OLTF’s  must  have  a  phase  margin  angle  of  at  least  30  degrees  and  a  gain  margin  of  at  least 
6  dB. 


5.2.3  Lateral  Performance  Specifications.  The  Level  1  flying  quality  roll  performance 
requirements  from  the  Mil  Std  are  used  as  the  lateral  FCS  performance  benchmarks.  These  bench¬ 
marks  must  be  met  in  order  to  achieve  an  FCS  comparable  to  the  current  F-16  FCS.  The  only 
difference  between  the  lateral  and  longitudinal  channel  benchmarks  is  that  the  lateral  benchmarks 
are  not  physically  attainable  by  aircraft  of  low  dynamic  pressure.  Therefore,  during  the  time  do¬ 
main  analysis,  the  dynamic  pressure  at  which  a  Level  1  roll  requirement  is  met  must  be  subjectively 
evaluated  as  to  its  acceptability.  The  roll  performance  benchmarks  are 

1.  Level  1  -  360°  of  bank  in  <  2.8  seconds 


5-7 


Figure  5.6  Lateral  3X2  MIMO  Plant  Structure  including  Dutch  Roll  Damping 


2.  Level  1  (combat)  -  90"  of  bank  in  <  1  second 


5.2.4  QFT  Lateral  Channel  Design. 


5.2.4. 1  QFTCAD  Setup.  Lateral  channel  MIMO  plants  used  in  the  lateral  QFT 
design  have  3  control  inputs  {Sdjtaih  ^aii,  and  and  two  outputs  (p  and  /?).  Matlab  is  used 

to  generate  the  appropriate  transfer  functions  to  represent  the  3X2  MIMO  system  illustrated  in 
Figure  5.6.  Note  that  the  previously  designed  dutch  roll  damping  loop  is  incorporated  into  the 
MIMO  system  plant  model. 

Once  generated,  the  282  lateral  plants  are  loaded  into  the  QFTCAD.  The  lateral  design 
structure  entered  into  the  QFTCAD  is  illustrated  in  Figure  5.7.  The  weighting  matrix  W  allocates 
Pcmd  into  and  and  transforms  the  3X2  MIMO  system  into  an  effective  2X2 

MIMO  system.  In  Reynold’s  research,  it  was  determined  that  the  weighting  matrix  of  Equation  (5.7) 
yielded  the  best  results  and  is  therefore  adopted  for  use  in  this  research  [8]. 


(5.7) 
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Figure  5.7  Lateral  Design  QFTCAD  Structure 


Next,  the  effective  plant  Pg  (basic  plant  plus  weighting  matrix  and  actuators)  and  q  transfer 
functions  are  generated  by  the  QFTCAD.  During  the  generation  of  these  transfer  functions,  the 
QFTCAD  cancels  any  poles  and  zeros  within  0.001  of  one  another  and  eliminates  all  poles  and  zeros 
outside  the  dynamic  frequency  range  of  10“®  to  10^  rps.  All  282  lateral  plant,  effective  plant,  and 
q  transfer  functions  are  listed  in  the  supplement  to  this  thesis.  Finally,  the  specifications  presented 
in  Section  5.2.2  are  entered  into  the  QFTCAD. 

Diagonal  Dominance.  The  diagonal  dominance  of  the  Q  matrices  is  evaluated 
by  the  QFTCAD  and  results  in  the  plots  of  Figure  5.8.  Diagonal  dominance  is  a  measure  of  the 
dominance  of  the  system’s  response  by  the  diagonal  MISO  loops  over  the  off-diagonal  MISO  loops. 
The  full  explanation  of  this  condition  is  found  elsewhere  [6].  At  the  frequencies  where  the  QFTCAD 
diagonal  dominance  plots  of  Figure  5.8  are  negative,  diagonal  dominance  does  not  exist;  therefore, 
diagonal  dominance  does  not  exist  in  the  design  frequency  spectrum  of  0.5  to  3.5  rps. 

Three  options  are  available  when  this  condition  of  diagonal  dominance  is  not  met.  The  first 
option  is  to  attempt  a  method  2  QFT  design.  In  a  method  2  design,  the  MISO  loop  expected  to 
have  the  largest  bandwidth  is  designed  first.  When  the  next  MISO  loop  is  designed,  the  cross¬ 
coupling  input  into  the  second  loop  reflects  the  compensator  design  of  the  first  loop.  In  contrast, 
when  using  method  1,  each  loop  is  designed  independently  assuming  the  worst  case  cross-coupling 
effect  possible  for  each  loop.  Therefore,  method  2  decreases  the  overdesign  inherent  in  the  second 


5-9 


Figure  5.8  QFTCAD  Diagonal  Dominance  Plot 


and  subsequent  MISO  loop  designs  that  occurs  when  using  method  1  [6] .  The  second  option  is  to 
modify  the  weighting  matrix  until  the  diagonal  dominance  condition  is  met.  Although  techniques 
for  arriving  at  such  a  weighting  matrix  exist,  they  are  not  useful  when  large  parameter  uncertainty 
is  present;  consequently,  trial  and  error  must  be  relied  upon  to  properly  modify  W.  The  third 
option  is  to  proceed  with  a  method  1  design,  but  to  design  with  only  stability  in  mind.  In  this  case, 
time  domain  simulation  must  be  relied  upon  to  verify  that  appropriate  tracking  and  cross-coupling 
performance  is  achieved. 

The  third  option  is  chosen  for  this  design  for  several  reasons.  The  lack  of  diagonal  dominance 
drastically  affects  the  generation  of  QFT  tracking  and  cross-coupling  rejection  boundaries.  In 
essence,  it  invalidates  these  boundaries.  A  method  2  design  is  attempted,  but  the  boundaries 
generated  by  the  QFTCAD  are  essentially  the  same  as  those  for  method  1.  In  other  words,  method 
2  does  not  rectify  the  diagonal  dominance  problem  and  the  first  option  is  a  dead  end.  Since  Reynolds 
attempted  to  modify  the  W  matrix  during  his  research  with  poor  results,  the  second  option  is  not 
pursued  [8].  A  final  consideration  favoring  the  third  option  is  that  the  detection  of  actuator  rate 
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Figure  5.9  QFT  Roll  Rate  Loop  Frequency  Templates 


and  control  surface  deflection  saturations  requires  time  simulation  of  the  system;  hence,  the  time 
simulations  can  easily  be  used  to  also  examine  the  qualitative  characteristics  of  the  system  response. 

5.2.4-S  QFT  Lateral  Roll  Rate  Loop  Design.  The  QFTCAD  frequency  templates 
for  the  roll  rate  loop  are  shown  in  Figure  5.9.  As  indicated  in  the  figure,  plant  number  276  has 
been  chosen  cis  the  nominal  plant  due  to  its  position  at  the  top  to  the  30  rps  template. 

The  roll  rate  nominal  loop  shape  is  depicted  on  the  Nichols  chart  of  Figure  5.10.  The  nominal 
loop  satisfies  all  of  the  stability  bounds  on  the  chart.  The  compensator  used  to  achieve  this  loop  is 


^  0.125(s  +  3) 

Gp  - 


(5.8) 


The  roll  rate  loop  prefilter  is  chosen  in  the  same  manner  as  that  for  the  longitudinal  control 
loop.  The  prefilter 


^  3.5 

"  “  s  +  3.5 


(5.9) 
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allows  pilot  roll  commands  of  up  to  3.5  rps  to  pass  to  the  compensator;  therefore,  the  pilot’s  control 
bandwidth  is  not  restricted  by  the  prefilter. 

5. 2. 4. 3  QFT  Lateral  P  Loop  Design.  The  QFTCAD  frequency  templates  for  the  P 
loop  are  shown  in  Figure  5.11.  As  indicated  in  the  figure,  plant  number  276  has  been  chosen  as  the 
nominal  plant  due  to  its  position  at  the  top  to  the  30  rps  template. 

The  P  nominal  loop  shape  is  depicted  on  the  Nichols  chart  of  Figure  5.12.  The  nominal  loop 
does  satisfy  all  of  the  stability  bounds  on  the  chart.  The  compensator  used  to  achieve  this  loop  is 


db 
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Figure  5.12  QFTCAD  ^  Nominal  Loop  Nichols  Chart 
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Plants  below  ( - ]  and  above  (solid)  150  psf  dynamic  pressure 


Figure  5.13  Initial  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (1  of  4) 

The  selection  of  a  /?  loop  prefilter  is  discussed  in  the  time  analysis  section  since  the  filtering  of 
pilot  sideslip  inputs  is  primarily  dependent  on  the  cross-coupling  roll  rate  observed  from  a  sideslip 
input. 


5.2.5  Lateral  Design  Time  Domain  Analysis.  Time  responses  to  a  unit  roll  rate  step 
input  are  made  at  this  point  as  a  progress  check  of  the  design.  The  Matlab  plots  resulting  from 

this  simulation  for  representative  plants  below  ( - )  and  above  (solid)  q  =  150  psf  are  shown  in 

Figures  5.13  through  5.16.  Two  deficiencies  of  the  design  are  revealed  by  inspecting  these  plots. 
First,  the  p  responses  of  the  low  q  plants  in  Figure  5.13  possess  an  undesirable  droop  prior  to 
reaching  the  desired  final  value.  Second,  as  seen  in  Figure  5.14,  the  low  q  plants  also  violate  the 
0.067  degree  /?  restriction  by  a  wide  margin. 

The  first  problem  is  rectified  by  decreasing  the  gain  of  the  r  feedback  washout  filter  of  Equa¬ 
tion  (5.1)  from  2.6  to  1.  The  rudder  deflection  commanded  by  the  washout  filter  with  a  gain  of  2.6 
slows  the  roll  rate  and  results  in  the  droop  in  the  p  response.  By  reducing  the  washout  filter  gain, 
thus  reducing  the  magnitude  of  the  rudder  deflection,  this  droop  is  adequately  reduced. 
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Figure  5.14  Initial  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (2  of  4) 
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Figure  5.15  Initial  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (3  of  4) 
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Plants  below  ( - )  and  above  (solid)  150  psf  dynamic  pressure 


Figure  5.16  Initial  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (4  of  4) 

The  second  problem  is  solved  by  increasing  the  gain  of  the  of  Equation  (5.10).  This  action 
reduces  the  sideslip  resulting  from  a  roll  rate  command  input  and,  because  /?  is  a  slower  variable 
than  r,  this  change  in  gain  has  much  less  effect  on  the  roll  rate  response  than  a  change  in  the  gain 
of  the  washout  filter. 

Modified  MIMO  plants  are  generated  by  Matlab  using  the  lower  washout  filter  gain  and  then 
input  into  the  QFTCAD.  The  Gp  of  Equation  (5.8)  is  unchanged,  and  the  absolute  gain  of  the 
Gp  of  Equation  (5.10)  can  now  be  increased  to  50  without  violating  stability  specifications.  The 
resulting  QFTCAD  nominal  OLTF’s  for  the  p  and  /?  loops  are  illustrated  in  Figures  5.17  and  5.18. 

Roll  rate  unit  step  input  time  responses  reflecting  the  adjusted  Gp  and  washout  filter  are 
depicted  in  Figures  5.19  through  5.22.  Figure  5.19  clearly  shows  that  the  p  response  droop  has  been 
reduced;  however,  the  maximum  transient  /?  of  the  low  q  plants  still  exceeds  the  0.067  degree  limit. 
Consequently,  the  0  compensator  design  must  be  divided  into  two  parts  because  the  gain  of  Gp 
cannot  be  increased  further  without  violating  the  gain  margin  restriction  of  6  dB  (see  Figure  5.18). 
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Plants  bslow  ( - )  and  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.19  Adjusted  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (1  of  4) 


Plants  below  ( - )  arxl  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.20  Adjusted  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (2  of  4) 
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Plants  below  ( — )  and  above  (solid)  ISO  psf  dynamic  pressure 
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Figure  5.21  Adjusted  Lateral  Design  Roll  Rate  Unit  Step  Input  Time  Response  (3  of  4) 


Plants  below  ( — }  and  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.23  QFTCAD  Expanded  Frequency  Template  for  ^  Loop 


5. 2. 5.1  Modified  fi  Loop  Design.  The  gain  of  Gp  needs  to  be  increased  for  the  low  q 
plants  in  order  to  comply  with  the  0.067  degree  sideslip  restriction.  First,  a  determination  is  made 
at  to  where  in  the  parameter  space  the  compensator  design  will  be  divided.  QFTCAD  expanded 
frequency  templates  such  as  the  one  depicted  in  Figure  5.23  are  examined  to  determine  a  prudent 
dynamic  pressure  at  which  to  divide  the  design.  Referencing  Figure  5.23,  the  division  is  made 
beween  plants  179  and  180  corresponding  to  g  =  150  psf. 

The  unit  p  step  input  time  response  of  the  low  q  plants  with  a  Gg  gain  of  -100  are  shown 
in  Figures  5.24  and  5.25.  Figure  5.25  validates  that  this  gain  increase  does  reduce  the  transient 
sideslip  to  less  than  0.067  degrees;  consequently,  the  gain  of  Gp  is  switched  between  -100  and  -50 
at  5  =  150  psf.  Below  a  g  of  150  psf,  a  gain  of  -100  is  used,  and  when  q  increases  above  150  psf,  the 
gain  of  Gp  is  changed  to  -50.  An  added  benefit  of  the  gain  switching  is  that  the  p  response  of  the 
low  q  plants  is  improved  over  that  shown  in  Figure  5.20.  The  p  response  ballooning  exhibited  by 
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Plants  below  150  psf  dynamic  pressure 


Figure  5.24  Unit  Roll  Rate  Step  Input  Response  of  Low  q  Plants  with 
(1  of  2) 


the  lowest  q  plant  of  Figure  5.24  is  isolated  to  plants  corresponding  to  dynamic  pressures  between 
60  and  70  psf. 

To  facilitate  the  required  gain  switching,  overlapping  of  the  dynamic  pressure  regions  accom¬ 
modated  by  the  two  compensator  gains  is  desired.  Therefore,  the  two  dynamic  pressure  regions 
of  interest  are  selected  as  q  <  170  psf  and  q  >  130  psf.  The  nominal  loops  for  the  low  and  high 
dynamic  pressure  regions  are  shown  in  Figures  5.26  and  5.27,  respectively.  As  seen  in  the  figures, 
both  nominal  loops  meet  the  phase  margin  angle  and  gain  margin  specifications. 


5. 2. 5. 2  Synthesis  of  Lateral  Maximum  Command  Input  Profiles.  The  roll  rate  loop 
is  evaluated  first  since  it  is  the  primary  control  loop  in  the  lateral  channel.  The  lateral  system 
is  subjected  to  a  unit  roll  rate  step  input  to  determine  if  sideslip  limitations  are  violated  and  to 
generate  bounds  on  Pcmd  to  achieve  desired  performance  while  avoiding  actuator  rate  and  control 
surface  deflection  saturations.  The  results  of  this  simulation  for  a  representative  set  of  plants 
spanning  the  entire  range  of  q  are  displayed  in  Figures  5.28  through  5.31.  A  complete  tabular 
listing  of  these  results  for  all  plants  is  found  in  Appendix  E.l.  Examination  of  the  figures  and  the 
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60. 


C^n  Loop  Trmntmi««ionc>(«)'*qii(«)  for 


Figure  5.27  Nominal  Open-Loop  Nichols  Chart  for  q  >  130  psf  and  Gp  =  ~ 

appendix  data  reveals  that  the  responses  of  all  plants  remain  below  the  maximum  /?  limit  of  ±0.067 
degrees. 

Next,  upper  and  lower  bounds  for  the  allowable  maximum  pemd  at  all  dynamic  pressures  are 
generated  from  the  unit  step  response  data  in  Appendix  E.l  by  plotting  the  following  versus  q. 


60 

(5.11) 

60 

(5.12) 

20 

(5.13) 

30 

(5.14) 

360 

^(2.8  sec) 

(5.15) 

90 

^(1  #cc) 

(5.16) 

5-23 


Plants  below  ( - )  and  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.30  Lateral  Design  Unit  Pcmd  Step  Responses  (3  of  4) 


Plants  below  ( - )  arxi  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.31  Lateral  Design  Unit  pcmd  Step  Responses  (4  of  4) 
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Figure  5.32  Proposed  Profile  to  Avoid  Saturations  and  Achieve  Desired  Performance 

The  plots  of  the  data  generated  for  all  plants  from  Equations  (5.11)  through  (5.16)  are  shown  in 
Figure  5.32  along  with  a  proposed  profile.  The  Pcmd(„„,)  profile  must  stay  below  the  plots 

of  Equations  (5.11)  through  (5.14)  to  avoid  actuator  rate  and  control  surface  deflection  saturations. 
In  addition,  the  profile  achieves  the  level  1  (yf^)  and  combat  level  1  {-{-—f  )  roll  performance 
requirements  when  it  is  above  the  bounds  defined  by  Equations  (5.15)  and  (5.16),  respectively.  The 
figure  indicates  that  the  proposed  profile  avoids  all  saturations,  achieves  level  1  roll  performance 
above  q  «  200  psf,  and  achieves  level  1  combat  roll  performance  above  q  «  300  psf. 

The  sideslip  channel  prefilter  is  now  synthesized  by  examining  the  cross-coupling  effects  re¬ 
flected  in  p  from  a  ^cmd  step  input.  The  time  responses  to  a  ^cmd  unit  step  are  displayed  in 
Figures  5.33  through  5.36  for  a  represetative  set  of  plants  using  a  j3  prefilter  of 


0.5 

s  -1-  0.5 


(5.17) 


This  prefilter  is  chosen  on  the  basis  that  fast  responses  to  /3  commands  from  the  pilot  are  not  a 
necessity  since  /?  is  rarely  commanded  by  the  pilot  to  accomplish  lateral  tracking  tasks.  Also,  the 
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Plants  below  ( - )  and  above  (solid)  150  psf  dynamic  pressure 


time  (sec) 


Figure  5.33  Lateral  Design  Unit  ficmd  Step  Responses  (1  of  4) 

frequency  responses  of  MISO  loop{2,l}  in  Figure  5.38  indicate  that  significant  cross-coupling  from 
a  input  to  roll  rate  are  present;  therefore,  the  ^  inputs  of  the  pilot  are  heavily  filtered  to  reduce 
the  effects  of  this  cross-coupling.  A  complete  tabulation  of  the  unit  ficmd  step  input  responses  for 
all  plants  is  found  in  Appendix  E.3. 

In  the  same  fashion  as  the  roll  rate  loop,  the  unit  0cmd  input  step  response  data  is  used  to 
synthesize  a  profile  that  avoids  rate  and  deflection  saturations.  To  facilitate  this,  the  data 

generated  for  all  plants  by  Equations  (5.11)  through  (5,14)  is  plotted  versus  dynamic  pressure  to 
form  an  upper  bound  for  the  profile.  As  long  as  the  profile  remains  below  the  composite 

bound  formed  by  the  minimum  of  Equations  (5.11)  through  (5.14)  for  all  q,  no  rate  or  deflection 
saturations  will  occur.  A  plot  of  the  composite  upper  bound  and  the  proposed  profile  is 

shown  in  Figure  5.37.  No  performance  lower  bound  is  applicable  since  no  performance  requirements 
are  dictated  for  /?  response. 

5.2.6  Lateral  Design  Validation. 
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Plants  below  ( - )  and  above  (solid)  150  psf  dynamic  pressure 
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Figure  5.34  Lateral  Design  Unit  Pcmd  Step  Responses  (2  of  4) 


Plants  below  ( - )  arxl  above  (solid)  150  psf  dynamic  pressure 


Figure  5.35  Lateral  Design  Unit  Pcmd  Step  Responses  (3  of  4) 
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PlantG  below  ( — )  and  above  (solid)  150  pet  dynamic  pressure 


Figure  5.36  Lateral  Design  Unit  pcmd  Step  Responses  (4  of  4) 
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Figure  5.37  Proposed  Profile  to  Avoid  Saturations 
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5.2.6. 1  QFTCAD  Frequency  Domain  Lateral  Design  Validation.  The  QFTCAD 
tracking  validation  closed-loop  Bode  plots  for  the  MIMO  system  are  displayed  in  Figures  5.38  and 
5.39.  The  four  loops  of  the  figures  correspond  to  the  following  input  output  relationships. 

•  MISO  loop{l,l}  -  output  from  0  input 

•  MISO  loop{l,2}  -  0  output  from  p  input 

•  MISO  Ioop{2,l}  -  P  output  from  0  input 

•  MISO  loop{2,2}  -  p  output  from  p  input 

Figure  5.38  depicts  the  tracking  response  for  the  low  q  plants.  Recall  that  the  tracking  bound¬ 
ary  models  displayed  for  MISO  loop{l,l}  are  arbitrarily  chosen  to  enable  the  input  of  a  /?  loop 
prefilter  into  the  QFTCAD.  The  synthesis  of  the  0  loop  prefilter  is  discussed  in  Section  5. 2. 5. 2. 
The  frequency  responses  of  MISO  loop{l,l}  are  robust,  but  do  display  some  cross-coupling  ef¬ 
fects  corresponding  to  the  peak  of  the  MISO  loop{l,2}  frequency  response  plots.  The  plots  for 
MISO  loop{l,2}  indicate  that  some  plants  violate  the  0.067  cross-coupling  bound,  so  compliance 
with  0  cross-coupling  limits  must  be  confirmed  through  time  domain  analysis.  MISO  loop{2,l}  of 
Figure  5.38  emerges  as  the  reason  for  the  lack  of  diagonal  dominance  in  the  lateral  channel.  Signif¬ 
icant  cross-coupling  of  sideslip  to  roll  rate  is  present,  especially  below  1  rps.  The  MISO  loop{2,2} 
frequency  responses  definitely  reflect  the  effect  of  the  cross-coupling,  but  are  deemed  acceptable. 
The  physical  explanation  for  the  siginificant  cross-coupling  is  that  the  dutch  roll  mode  and  0  regu¬ 
lation  loops  impact  the  roll  rate  response  at  low  dynamic  pressures.  As  a  pilot  initiates  a  roll  at  low 
dynamic  pressures,  the  rudder  input  commanded  by  the  dutch  roll  mode  damping  and  0  regulation 
loops  slows  the  roll  rate.  This  results  in  the  attenuation  and  delay  of  the  roll  rate  response  in  the 
frequency  domain.  The  closed  loop  bandwidth  of  all  of  the  low  q  plants  is  less  than  3.5  rps,  but 
this  fact  makes  physical  sense  since  the  aircraft  is  less  responsive  at  these  lower  energy  states. 

In  contrast  to  the  responses  of  the  low  q  plants,  the  high  q  plant  responses  illustrated  in 
Figure  5.39  are  very  robust  and  MISO  loops{l,l}  and  {2,2}  show  only  slight  evidence  of  cross- 
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MISO  loop  {i.l) 


MISO  loop{l^} 


coupling  effects.  All  of  the  responses  of  MISO  loop{l,2}  remain  below  the  desired  cross-coupling 
bound,  and  the  responses  of  MISO  loop{2,l}  show  that  the  cross-coupling  effects  on  roll  rate  are 
much  reduced  at  higher  dynamic  pressures. 


The  stability  validations  for  the  roll  rate  and  the  high  and  low  q  sideslip  loops  are  depicted 
on  the  Nichols  charts  in  Figures  5.40,  5.41,  and  5.42,  respectively.  These  figures  indicate  that  the 
required  stability  has  been  achieved  by  both  lateral  loop  designs. 


5.2.6.S  Lateral  Design  Time  Domain  Validation.  The  step  input  time  responses  of 
the  lateral  system  to  the  profile  are  depicted  in  Figures  5.43  through  5.46  for  a  repre¬ 

sentative  set  of  plants.  A  table  containing  the  results  of  this  simulation  for  all  plants  is  contained 
in  Appendix  E.2.  Figures  5.43  through  5.46  along  with  the  appendix  data  clearly  show  that  the 
sideslip  limit  of  ±6  degrees  is  not  exceeded  and  that  no  actuator  rate  or  control  surface  deflection 
saturations  are  present. 
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MISO  loop  (1.1} 


MISO  loop  (1.2) 


Figure  5.40  QFTCAD  Stability  Validation  for  Roll  Rate  Loop 
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Figure  5.41  QFTCAD  Stability  Validation  for  High  q  Sideslip  Loop 


Figure  5.42  QFTCAD  Stability  Validation  for  Low  q  Sideslip  Loop 
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Figure  5.44  Lateral  Design  Step  Responses  (2  of  4) 


Plants  below  ( — )  and  above  (solid)  1 50  pel  dynamic  pressure 


time  (sec) 


Figure  5.45  Lateral  Design  Pcmd(„.j)  Step  Responses  (3  of  4) 


Plants  below  (— )  and  above  (solid)  1 50  pel  dynamic  pressure 
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Figure  5.46  Lateral  Design  Pcmd(„„)  Step  Responses  (4  of  4) 
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Plants  balow  ( - )  and  above  (solid)  150  psf  dynamic  pressure 


time  (sec) 


Figure  5.47  Lateral  Design  Step  Responses  (1  of  4) 

The  step  input  time  responses  of  the  lateral  system  to  the  profile  are  depicted  in 

Figures  5.47  through  5.50  for  a  representative  set  of  plants.  A  table  containing  the  results  of  this 
simulation  for  all  plants  is  contained  in  Appendix  E.4.  These  figures  and  the  appendix  data  confirm 
that  no  actuator  rate  or  control  surface  deflection  limits  are  exceeded. 
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Plants  below  ( — )  and  above  (solid)  150  pst  dynamic  pressure 
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Figure  5.48  Lateral  Design  Step  Responses  (2  of  4) 


Plants  below  ( - and  above  (solid)  1 50  psf  dynamic  pressure 


time  (sec) 


Figure  5.49  Lateral  Design  Step  Responses  (3  of  4) 
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Plants  below  ( - )  and  above  (solid)  150  psi  dynamic  pressure 


time  (sec) 

Figure  5.50  Lateral  Design  Step  Responses  (4  of  4) 


5.3  Final  Lateral  Channel  Design 

The  final  design  in  the  lateral  channel  is  best  described  by  the  block  diagram  illustrated  in 
Figure  5.51.  To  summarize,  the  roll  stick  force  from  the  pilot  is  first  scaled  by  the  maximum  pcmd 
profile  and  then  filtered.  The  filtered  input  is  then  passed  to  the  roll  rate  compensator  feedback 
loop  where  it  is  transformed  into  an  appropriate  aileron  command.  The  horizontal  tails  receive  a 
differential  command  equal  to  approximately  one-third  of  the  aileron  command.  A  rudder  force 
input  from  the  pilot  is  first  scaled  by  the  maximum  Pcmd  profile  and  then  filtered.  The  filtered 
input  is  then  passed  to  the  0  compensator  feedback  loop  to  be  transformed  to  a  rudder  reference 
signal.  This  rudder  reference  signal  is  then  adjusted  by  the  dutch  roll  mode  damping  feedback 
loop  to  generate  an  appropriate  rudder  command.  The  roll  rate  prefilter  and  both  compensators 
meet  the  design  guidelines  set  forth  in  Section  2.3.  The  /3  prefilter  does  not  comply  with  the  3.5 
rps  minimum  bandwidth  requirement  for  the  reasons  discussed  in  Section  5. 2. 5. 2.  The  necessary  /? 
channel  compensator  gain  switching  at  g  =  150  psf  is  clearly  indicated  in  Figure  5.51.  This  is  the 
only  direct  gain  scheduling  present  in  either  the  longitudinal  or  lateral  channel  designs. 
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Figure  5.51  Final  Lateral  FCS  Design  Block  Diagram 
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VI.  Conclusions  and  Recommendations 

6.1  Conclusions 

In  line  with  the  research  objectives  proposed  for  this  thesis,  a  full  subsonic  envelope  FCS  has 
been  designed  for  the  VISTA  F-16.  An  integrated  approach  utilizing  QFT  along  with  time  domain 
simulation  has  resulted  in  a  design  that  is  not  only  quite  simple,  but  one  that  clearly  reflects  the 
physical  limitations  imposed  by  the  “real  world”  (i.e.,  acutator  rates  and  control  surface  deflections). 
QFT  was  an  invaluable  tool  in  the  practical  and  eflScient  achievement  of  this  task,  allowing  the 
designer  to  make  insightful  modifications  to  the  design  with  full  knowledge  of  the  consequences  of 
these  modifications. 

The  importance  of  staying  as  close  to  the  “real  world”  as  possible  when  utilizing  any  robust 
control  design  technique  cannot  be  overstated.  In  both  the  longitudinal  and  lateral  channels, 
the  first-cut  QFT  compensator  designs  produced  very  robust  time  responses  to  small  command 
inputs;  however,  when  inputs  of  the  magnitude  required  to  achieve  the  desired  performance  were 
considered,  lengthy  actuator  rate  and  control  surface  deflection  saturations  occurred.  Therefore, 
the  robustness  achieved  by  the  high  gain  of  these  initial  compensators  had  to  be  exchanged  for  the 
desired  performance  requirements  by  reducing  the  compensator  gain.  As  an  illustration,  Reynold’s 
final  body  axis  roll  rate  compensator 

_  1.6(s-b7.6)(s-H0) 

^  s{s  +  100) 


(6.1) 


of  this  thesis  which  has  a  DC  gain  of  0.375  [8].  The  responses  of  Reynold’s  system  were  very 
robust  in  response  to  small  command  inputs;  however,  with  a  compensator  gain  more  than  three 
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times  that  necessary  to  meet  design  tracking  requirements  while  avoiding  saturations,  Reynold’s 
suspicion  about  the  possiblity  of  saturations  when  the  system  is  acted  upon  by  large  inputs  is 
clearly  confirmed  [8]. 

This  FCS  is  void  of  any  direct  compensator  gain  scheduling  except  in  the  sideslip  control 
loop.  However,  indirect  scheduling  of  a  sort  is  incorporated  in  all  control  loops  through  the  scaling 
of  the  control  stick  force  inputs  of  the  pilot  by  a  g  dependent  profile.  These  maximim  allowable 
command  input  profiles  are  the  means  by  which  actuator  rate  and  control  deflection  saturations  are 
averted  in  this  design,  while  at  the  same  time  achieving  the  desired  performance.  In  this  framework 
(and  with  a  specified  prefilter),  the  compensator  gains  themselves  serve  to  establish  the  size  of  the 
allowable  range  of  the  maximum  command  input  profiles.  If  the  compensator  gain  is  too  high,  then 
the  actuator  rate  and  control  deflection  upper  bounds  fall  toward  the  performance  bounds,  thus 
closing  the  window  for  an  acceptable  maximum  command  input  profile.  An  example  of  this  effect 
of  compensator  gain  is  illustrated  in  Figure  6.1.  Although  higher  gain  increases  robustness,  it  could 
result  in  lower  than  desired  performance  once  saturations  are  taken  into  account.  The  location  of 
the  poles  and  zeros  of  the  compensators  affect  the  qualitative  shape  of  the  responses  while  also 
ensuring  stability. 

The  application  of  the  specifications  set  forth  in  Military  Standard  1797A  to  an  FCS  design 
proves  to  be  a  difficult  task.  Many  of  the  specifications  are  open-ended  or  incomplete.  Some 
specifications  are  defined  in  both  the  time  and  frequency  domains,  but  have  no  evident  relationship 
to  one  another.  Others  are  defined  by  only  an  upper  or  lower  limit,  rather  than  by  both  an  upper 
and  lower  boundary  of  acceptable  performance.  Consequently,  the  engineer  is  given  no  clear-cut 
choice  of  what  specifications  to  apply  to  a  design.  The  validity  of  the  particular  specifications  that 
are  chosen  to  guide  this  design  can  only  be  tested  by  the  evaluation  of  the  resulting  flying  qualities 
by  an  actual  pilot. 
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Upper  Saturation  Bound  ( - ) 


Figure  6.1  Effect  of  Compensator  Gain  Increase  on  Command  Input  Profile  Bounds 


The  success  of  imbedding  flying  qualities  into  the  longitudinal  design  by  utilizing  the  dynamic 
C*  control  parameter  can  only  be  measured  by  simulating  the  FCS  with  a  pilot  in  the  loop. 
Although  the  blend  of  plant  parameters  used  in  this  thesis  to  create  C*  may  not  be  optimal,  it 
provides  an  adequate  starting  point  from  which  to  continue  the  development  of  the  concept. 


The  lack  of  diagonal  dominance  in  the  lateral  channel  resulted  in  applying  QFT  to  only  the 
stability  specifications,  leaving  the  tracking  specifications  up  to  time  domain  simulation  analysis. 
This  resulted  in  more  trial  and  error  than  was  necessary  in  the  longitudinal  design,  but  the  QFT- 
CAD  was  still  instrumental  in  allowing  quick  adjustments  to  the  compensator  to  achieve  the  desired 
tracking  performance  while  complying  with  stability  constraints. 


The  FCS  of  this  research  effort  has  been  validated  thoroughly  in  both  the  time  and  frequency 
domains.  In  the  lateral  and  longitudinal  channels,  some  low  q  plants  exhibit  ballooning  in  response 
to  roll  rate  and  C*  command  inputs.  Whether  or  not  this  ballooning  constitutes  a  nuisance  to  a 
pilot  can  only  be  determined  by  simulating  the  FCS  with  an  actual  pilot  in  the  loop.  However,  if 
the  ballooning  is  deemed  a  problem,  it  can  be  rectified  by  splitting  the  roll  rate  and  C*  compensator 
designs  as  is  done  with  the  sideslip  loop  to  achieve  more  robustness  at  low  dynamic  pressures.  This 
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action  will  result  in  the  degredation  of  large  signal  input  performance  due  to  the  higher  compensator 
gain  needed  to  achieve  the  robustness. 

The  time  domain  validation  of  this  design  is  based  on  step  input  responses.  Other  inputs  that 
may  be  of  interest  in  FCS  design  include  pulses  and  doublets.  Pulse  or  doublet  control  inputs  may 
cause  higher  transient  actuator  rates  and  control  surface  deflections  than  those  of  a  step.  If  this  is 
the  case,  the  gain  of  the  compensator  and/or  the  maximum  command  input  profile  can  be  adjusted 
to  eliminate  the  saturations.  The  end  result  will  be  a  slight  degredation  of  aircraft  performance  in 
exchange  for  the  elimination  of  any  destabilizing  saturations. 

To  summarize,  this  research  effort 

1.  Introduces  a  new  concept  of  imbedding  flying  qualities  in  an  FCS  through  the  use  of  a  dynamic 
control  parameter  which  attempts  to  accurately  capture  the  true  desires  of  a  pilot  as  aircraft 
flight  conditions  vary. 

2.  Illuminates  the  benefits  of  applying  QFT  to  a  problem  of  large,  but  structured  uncertainty 
(full  subsonic  envelope  FCS  design). 

3.  Demonstrates  that  robust  control  design  techniques  can  be  applied  to  practical  problems 
and  result  in  physically  realizable  solutions,  but  only  when  limitations  imposed  by  the  “real 
world”  are  kept  close  at  hand  during  the  design  process. 

6.2  Recommendations 

To  fully  validate  the  concept  of  imbedding  flying  qualities  in  a  longitudinal  FCS  design  through 
the  use  of  the  C*  control  parameter  used  in  this  thesis,  it  is  recommended  that  a  full  nonlinear  sim¬ 
ulation  of  the  FCS  be  conducted.  If  the  nonlinear  simulation  does  not  reveal  significant  problems, 
the  FCS  should  then  be  simulated  with  a  pilot  in  the  loop.  The  flying  quality  assessment  resulting 
from  the  pilot-in-the-loop  simulation  will  determine  the  overall  effectiveness  of  the  C*  parameter 
and  will  also  allow  qualitative  assessment  of  the  flight  control  system  as  a  whole. 
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The  integrated  QFT  and  time  domain  simulation  method  utilized  in  this  design  required 
iteration  to  produce  a  viable  design.  If  a  way  of  generating  Nichols  chart  boundaries  to  reflect 
actuator  rate  and  control  surface  deflection  limitations  existed,  then  QFT  could  be  applied  without 
a  requirement  for  time  simulation.  These  physical  restriction  boundaries  might  appear  as  upper 
bounds  on  the  Nichols  chart.  In  this  case,  the  job  of  the  designer  would  be  to  synthesize  a  com¬ 
pensator  to  shape  the  nominal  open-loop  so  that  the  nominal  loop’s  magnitude-phase  points  are 
above  the  tracking  bounds,  but  below  these  physical  restriction  bounds  for  all  frequencies  in  the 
frequency  spectrum  of  the  design.  If  such  a  method  could  be  incorporated  into  the  QFTCAD,  the 
iteration  required  to  achieve  a  valid  design  could  be  greatly  reduced.  Therefore,  it  is  recommended 
that  future  research  efforts  explore  this  possibility. 

While  on  the  subject  of  pilot-in-the-loop  simulation,  further  research  directed  towards  incor¬ 
porating  pilot  compensation  into  the  FCS  design  process  is  warranted.  The  imbedding  of  flying 
qualities  in  control  parameters  or  placing  pilot  characteristics  in  system  prefilters  are  but  two  ex¬ 
amples  of  possible  types  of  pilot  compensation.  Such  research  cannot  hope  to  eliminate  the  need 
for  pilot-in-the-loop  simulation,  but  can  reduce  the  number  of  design  iterations  required.  In  other 
words,  if  an  FCS  is  designed  with  not  only  performance  specifications  but  also  pilot  characteristics 
and  desires  in  mind,  the  resulting  design  should  require  only  slight  modification  following  pilot-in- 
the-loop  simulation.  Therefore,  if  an  effective  means  of  emulating  the  characteristics  and  desires 
of  a  pilot  can  be  employed  from  the  outset  of  an  FCS  design,  the  design  process  will  certainly  be 
more  efficient  and  less  costly. 

The  Mil  Std  and  other  documents  like  it  that  are  used  to  guide  the  design  of  flight  control 
systems  need  to  continue  to  evolve  towards  the  definition  of  specifications  that  can  be  easily  in¬ 
corporated  at  the  onset  of  a  design  process  rather  than  those  that  must  be  be  evaluated  at  the 
completion  of  the  design.  The  task  of  developing  one  set  of  specifications  for  the  full  envelope  of 
an  aircraft  is  an  impossible  task  due  to  the  large  variations  in  aircraft  control  response;  however. 
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specifications  in  the  form  of  upper  and  lower  boundaries  of  acceptable  performance  could  certainly 
be  developed  for  specific  regions  of  the  flight  envelope. 

A  couple  of  suggested  improvements  to  the  QFTCAD  are 

1.  The  QFTCAD  currently  has  no  feature  to  facilitate  the  input  of  a  large  number  of  plant  mod¬ 
els.  This  feature  would  have  been  a  blessing  in  this  design  effort,  and  should  be  incorporated 
in  future  versions. 

2.  In  some  cases,  the  designer  is  interested  in  satisfying  QFT  tracking  or  disturbance  boundaries 
for  only  some  of  the  selected  template  frequencies.  In  fact,  the  frequencies  of  interest  for 
tracking  may  be  different  than  those  for  disturbance  rejection.  Currently,  the  QFTCAD 
composite  boundaries  are  constructed  using  the  tracking  and  disturbance  boundaries  for  all 
template  frequencies.  It  is  suggested  that  the  composite  boundary  feature  be  modified  to  allow 
the  selection  of  which  bounds  are  to  be  used  for  the  composite  at  each  template  frequency. 
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Appendix  A.  Augmented  Plant  Generation  and  QFTCAD  Plant  Input 

All  of  the  augmented  plants  used  in  this  thesis  are  generated  in  the  following  manner  to  avoid 
having  to  input  them  individually  into  the  QFTCAD. 

1.  The  augmented  system’s  block  diagram  is  constructed  in  the  Sima/ini  environment  [10]. 

2.  The  Matlab  ‘linmod’  command  is  used  to  generate  state  space  matricies  (A,  B,  C,  and  D)  from 
the  appropriate  input(s)  to  the  appropriate  output(s)  of  the  augmented  system  [11].  These 
matricies  are  generated  for  each  of  the  basic  SRF  longitudinal  or  lateral  aircraft  models. 

3.  The  state  space  matricies  are  then  factored  into  gain-zero-pole  form  using  the  ‘ss2zp’  Matlab 
function  [11].  This  is  the  form  that  the  QFTCAD  requires. 

4.  The  gain-zero-pole  transfer  functions  are  written  to  a  file  in  the  data  format  of  the  QFTCAD 
‘qftsave.m’  file  using  the  Matlab  ‘fprintf’  command  [11]. 

5.  The  formatted  file  is  pasted  into  a  ‘qftsave.m’  file  which  can  then  be  loaded  by  the  QFTCAD 
‘load’  feature  [1]. 

Concern  was  raised  as  to  the  validity  of  generating  the  augmented  plants  in  this  fashion  using 
the  Matlab  ‘linmod’  function.  The  plant  generation  was  validated  by  simulating  both  the  original 
system  block  diagram  and  one  containing  only  the  augmented  state  space  matricies.  The  time 
responses  of  both  the  original  and  augmented  system  representations  were  identical;  hence,  this 
plant  generation  technique  is  deemed  valid. 

An  example  Matlab  m-file  which  automates  this  entire  process  for  a  3  X  2  MIMO  control 
system  with  no  external  disturbance  follows  [11]. 
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%  Open  transfer  function  data  file 


fid  =  fopen(’filename’,’a’); 
fprintf(fid, ’plants  =  {’); 

%  Generate  gain-zero-pole  transfer  functions  from  the  SRF  state  space  plants 

for  u  =  1:282; 

eval([’m’  num2str(u)]); 

A=A.VISTA(5:8,5:8); 

B=[-B.VISTA(5:8,2),-B-VISTA(5:8,4),-B.VISTA(5:8,5)]; 

[a,b,c,d]=linmod(’5ima/inl;  model  name’); 

[zl,pl,kl]=ss2zp(a,b,c,d,l); 

[z2,p2,k2]=ss2zp(a,b,c,d,2); 

[z3,p3,k3]=ss2zp(a,b,c,d,3); 

zll=:[zl(:,l)]; 

z21=[zl(:,2)]; 

zl2=[z2(:,l)]; 

z22=[z2(:,2)]; 

zl3=[z3(:,l)]; 

z23=[z30,2)]; 

pll=[pl]; 

p21=[pl]; 

pl2=[p2]; 

p22=[p2]; 

Pl3=[p3]; 

p23=[p3]; 

kll=kl(l,l); 

k21=kl(2,l); 

kl2=k2(l,l); 

k22=k2(2,l); 

kl3=:k3(l,l); 

k23=k3(2,l); 

%  Write  the  gain,  zeros,  and  poles  to  the  data  file 

fid  =  fopen(’filename’,’a’); 
fprintf(fid,’{“%i”,  {},  {{’,u); 

%  First  output  from  first  input 

fprintf(fid,’transferF[s,  s,  0][%26.18E,\n  {’,kll); 
w=size(zll,l); 

if  fix(w/2)<(w/2); 

for  x=l:fix(w/2); 

fprintf(fid,’%26.18f-f%26.18PI,%26.18f-%26.18f*=I,\n’,real(zll(2*x-l)), 

abs(imag(zll(2’*'x-l))),real(zll(2*x)),abs(imag(zll(2*x)))); 

end; 

fprintf(fid,’%26.18f-|-%26.18f*’I},  {},  {’,real(zll(w)),abs(imag(zll(w)))); 
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for  x=l:(fbc(w/2)-l); 

fprintf(fid,’%26.18f+%26.18f*I,%26.18f-%26.18f*I,\n’,real(zll(2*x-l)), 

abs(imag(zll(2*x-l))),real(zll(2*x)),abs(imag(zll(2*x)))); 

end; 

fprintf(fid,’%26.18f+%26.18f*I,%26.18f-%26.18f*I},{},  {’,real(zll(w-l)), 
abs(imag(zll(w-l))),real(zll(w)),abs(imag(zll(w)))); 

end; 

w=size(pll,l); 

if  fix(w/2)<(w/2); 

for  x=l:fix(w/2); 

fprintf(fid,’%26.18f+%26.18PI,%26.18f-%26.18f^I,\n’,real(pll(2*x-l)), 

abs(imag(pll(2*x-l))),real(pll(2*x)),abs(imag(pll(2*x)))); 

end; 

fprintf(fid,’%26.18f+%26.18f*^I}.  {}],\n’,real(pll(w)),abs(imag(pll(w)))); 

else 

for  x=:l:(fix(w/2)-l); 

fprintf(fid,’%26.18f+%26.18PI,%26.18f-%26.18f*'I,\n’,real(pll(2*x-l)), 

abs(imag(pll(2*x-l))),real(pll(2*x)),abs(imag(pll(2*x)))); 

end; 

fprintf(fid,’%26.18f+%26.18Pl,%26.18f-%26.18m},{}],\n’,real(pll(w-l)), 

abs(iniag(pll(w-l))),real(pll(w)),abs(iniag(pll(w)))); 

end; 

%  First  output  from  second  input 

Same  as  for  ‘First  output  from  first  input’  with  “zll”,  “pH”,  and  “Atll”  replaced  by  “zl2”,  “pl2”, 
and  “ifcl2”. 


%  First  output  from  third  input 

Same  as  for  ‘First  output  from  first  input’  with  “zll”,  “pll”,  and  “fcll”  replaced  by  “zl3”,  “pl3”, 
and  “/fcl3”.  End  with  {}]},  {”  instead  of  {}],”. 

%  Second  output  from  first  second  and  third  input 

Same  as  for  first  output  with  all  “zl”,  “pi”,  and  “fcl”  replaced  by  “z2”,  “p2”,  and  “1:2”.  End 
with  {}]}},  Null},”  instead  of  {}]}, 

end; 

fclose(fid); 

%  Final  editing  of  ‘filename’  before  pasting  into  ‘qftsave.m’. 

%  A  UNIX  editor  source  file  is  created  which  changes  the  last  line  of  the  file  from  “...Null},” 
%  to  “...Null}}”  and  replaces  all  exponential  operators  “E”  with  “*10*”. 

fid=fopen(  ’edit  ’ ,  ’a’) ; 

fprintf(fid,’ed  <<  !\ne  filename\n,s/E/*10''/\n\\$s/Null},/Null}}/\nw  filename\nq\n!\n’); 
fclose(fid); 

Isource  edit 
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If  the  particular  system  of  interest  does  have  an  external  disturbance,  then  the  transfer  func¬ 
tions  from  the  disturbance  input  to  the  system  outputs  are  written  in  the  ‘filename’  file  in  place 
of  “Null”. 
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Appendix  B.  Parameter  Space  Data  Points 

This  appendix  contains  a  list  of  all  of  the  data  points  for  which  LTI  aircraft  models  are 
generated  by  the  SRF. 
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Appendix  C.  Longitudinal  Channel  Time  Domain  Simulation  Data 

This  appendix  contains  all  longitudinal  channel  time  domain  simulation  results. 
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C.l  Longitudinal  Unit  Step  Input  Time  Response  Data 
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185  699.9533  0.090  0.013  0.165  0.643  0.03137  5.010  0.453  1.765  1.195 

186  699.9533  0.090  0.013  0.165  0.643  0.03803  5-006  0.453  1.785  1.194 

187  771.1717  0.085  0.010  0.135  0.533  0.00076  4.615  0.396  1.838  0.964 

188  771.1717  0.085  0.010  0.135  0.533  0.00407  4.633  0.399  1.838  0.936 

189  771.1717  0.085  0.010  0.135  0.533  0.00311  4.836  0.398  1.838  0.936 

190  835.3739  0.060  0.009  0.130  0.516  0.00030  4.761  0.381  1.847  0.639 

191  635.3739  0.060  0.009  0.130  0.516  0.00006  4.758  0.360  1.847  0.839 

193  954.4603  0.080  0.009  0.130  0.516  0.00000  4.644  0.353  1.846  0.793 

193  954.4603  0.060  0.009  0.130  0.516  0.00031  4.664  0.356  1.846  0.769 

194  954.4603  0.080  0.009  0.130  0.516  0.00006  4.660  0.356  1.846  0.766 

195  1106.3130  0.075  0.009  0.115  0.511  0.00000  4.535  0.337  1.850  0.700 

196  1106.3130  0.075  0.009  0.115  0.511  0.00000  4.546  0.330  1.850  0.680 

197  1106.3130  0.075  0.009  0.115  0.511  0.00000  4.543  0.339  1.850  0.679 

198  1157.0650  0.075  0.009  0.110  0.499  0.00000  4.493  0.319  1.856  0.633 

199  1157.0650  0.075  0.009  0.110  0.499  0.00000  4.489  0.319  1.856  0.633 


C.2  Longitudinal  Maximum  C*cmd  Step  Input  Time  Response  Data  Without  Limiters 
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Appendix  D.  Rate  Limiter  Two-Dimensional  Arrays 


The  two-dimensional  array  used  to  limit  a  indexed  on  a  and  d  is, 


The  two-dimensional  array  used  to  limit  gpu  indexed  on  gpn  and  gpn  is, 


D-1 


Appendix  E.  Lateral  Channel  Time  Domain  Simulation  Data 
This  appendix  contains  all  lateral  channel  time  domain  simulation  results. 
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Time  Response  Data 
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E.4  Lateral  Maximum  Pcmd  Step  Input  Time  Response  Data 
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